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Abstract 

We present a summary of the basic properties of the radio wave generation, prop- 
agation and reception, with a special attention to the gigahertz bandwidth region 
which is of interest for wireless sensor networks. 



1 Introduction 



Over a period of several months we have made measurements with a set 
of transceivers with the purpose of investigating how the received power 
varies with the surrounding and placement of the devices. The RF-devices 
automatically measure a parameter called RSSI for Received Signal Strength 
Indicator, and thus provide a convenient means to track the power level of 
the signal. Since our measurements raised many issues about the behaviour 
of electromagnetic fields, it was decided to review some of the basic physics of 
electromagnetism in the style of a handbook chapter. The emphasis here is 
on the application of the Maxwell equations to concrete problems, not on the 
development of the theoretical structure (in terms of differential forms, gauge 



theory, etc). Of physics books on EM theory we may mention [19|, |26j, |36j 
and the engineering style books 42, 47], 13] • A good general physics reference 
including material on EM is [2lJ|. For applications to antenna theory see 28]. 
General reviews of EM with wireless networks in mind can be found in [l^E]. 
Arnold Sommerfeld was a an eminent mathematical physicist who made, 
among other things, some significant contributions to the propagation of EM 



fields; these "classical" methods are described in 14, 44 



45, 46 . Of recent 



texts on antenna theory we may mention (1, [3]. Balanis has also written a 
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nice minireview |2|. For an interesting online collection of lecture notes and 
a selection of classic papers on EM see [29 . 



2 Maxwell equations 

2.1 Quasi-stationary fields 

Electromagnetism stands for one of the four fundamental forces in physics. A 
static point like charge q\ at the point ri in an isotropic homogeneous medium 
exerts a force on an other charge g 2 at r 2 given by (Coulomb interaction), 



Fi_ 



1 gig 2 (r 2 -ri) 



47ree |r 2 - r : 

The quantity e, relative permittivity, is a 
quantity characterizing the medium, while eo 
(the permittivity of the vacuum) is a universal 
constant. Eq.fpQ) can be written as 



(1) 



Fi_* 2 = g 2 E(r 2 ), 



where 



E(r 2 ) 



1 gi(r 2 -ri) 
47ree |r 2 - ri| 3 



(2) 



(3) 




E(r 2 



is defined as the electric field at the point 
r 2 generated by the charge q\ located at i^. 

The electric field can also be expressed in terms of a potential function 



E 



-V0, 



or vice versa, 



0(r) = - / E-dr, (4) 

where the line integral is along a path connecting the reference point r 
and the point r. In the case of Eq.((3|) we have 

1 Qt 



Aireeo |r 2 - r x | 

When charges are in a relative motion with respect to each other then we 
have to include in Eq.([2|) a term depending on the velocity, 
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F = qE + gv x B, (5) 

where B defines the magnetic field strength. Thus, Eq.((5]) gives the force 
("Lorentz force") acting on a charge q moving with velocity v in an EM field 
characterized by E and B. From this follows the familiar fact, that a straight 
conductor of length I, with the current I in a magnetic field B, will sense a 
force BIl in a direction perpendicular to B and the conductor. 

Conversely, a charge q\ at r 1 moving with the velocity v x generates a 
magnetic field strength at the point r 2 given by (in an isotropic homogeneous 
medium) 

B(r 2 ) = ^ flvi*fo-r0 . (6) 
4tt |r 2 - ri| 3 v ' 

The field strength is thus affected by the magnetic permeability jj, char- 
acterizing the medium. For most non-metallics \i ~ 1, while /io is universal 
constant. From the above it follows that two charges moving with velocities 
vi and v 2 will interact via a magnetic force given b)Q 



T7i -p / \ gig2/x /x v 2 x (vi x (r 2 - n)) 

Fi^ 2 = g 2 v 2 x B(r 2 ) = — , R = (7) 

47r |r 2 — 

qiq 2 ^ (v 2 • (r 2 - ri))vi - (vi ■ v 2 )(r 2 - ri) 



Air |r 2 — ri| 3 

where we have used the rule that 



A x (B x C) = (A • C)B - (A ■ B)C. 

A current J in a conductor consists of many moving charges, each one 
contributing to the total magnetic field strength according to ((7j) . If we 
consider a small segment d\ of the conductor, then the sum of all terms qv 
over the charges in this segment is equal to / d\. Therefore the magnetic field 
strength generated by this segment is given by ("Biot-Savart law") 

1 An interesting observation is that the force Fi^ 2 which the particle 1 exerts on the 
particle 2 is no longer, in general, the opposite of the force that the particle 2 exerts on 1, 
as would be demanded by the principle of actio est reactio of Newtonian mechanics; that 
is, we no longer have Fi^ 2 + F 2 ^i = 0. From Fi^ 2 + F 2 ^i 7^ one might conclude that 
the closed system of charge 1 + charge 2 may start to move without an external cause. 
However, the momentum of the total system is conserved if we also take into account 
the momentum contribution of the electromagnetic field. Calculating the magnetic forces 
between two current loops, on the other hand, we get Fi_^ 2 + F 2 ^i = 0, in accordance 
with Newton's third law. 
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dB(r 2 ) 



Hofildl x (r 2 - ri) 



47r |r 2 — ri| 

that is, moving charges forming a current I 
in a small conducting element d\ at generates 
a magnetic field dB(r 2 ) given by Eq.® at the 
point r 2 . In order to obtain the effect of the 
whole conductor one has to sum (integrate) (El) 
over all the segments. 

The permittivity e and permeability take 
into account how the medium affects the elec- 
tromagnetic field. The charges in the medium 
are affected by the field and may become dis- 
placed, which leads to a modification of the 
field ("backreaction"). This explains such phe- 
nomena as polarization (charge displacements) 
and magnetization of a medium. From Eq.([3]) 
and Eq.([8]) we infer that by defining 



(8) 



I d\ 




rfB(r 2 



D = ee E, (9) 
H = — B, 

we obtain the quantities D (electric displacement) and H (magnetic field), 
which are apparently independent of the material factors (e, /i). From the 
definitions one can show that the integral of D over a boundary dV enclosing 
a volume V is equal to the total charge Q contained in V, while integrating 
H along a loop (boundary) dS enclosing a surface S one obtains the total 
current / flowing through that surface, 

/ B-dS = Q, (10) 

JdV 

I H-ds = I. 

JdS 

2.2 General case — time dependent fields 

If we integrate B over a surface S we obtain a quantity 

$ = f B • dS 

Js 



2 Maxwell equations 



6 



termed the magnetic flux through the surface S. It is experimentally 
observed that when the flux enclosed by a conducting loop changes, this 
induces a potential difference along the loop and causes a current to flow. 
More precisely (Faraday's law if induction), 



~dt = 



s dt 



E-ds. 



c)S 



Using the mathematical identity (Stokes' theorem) 



A • ds 



VxA-dS, 



'as js 
we obtain the induction law on the form, 



V x E 



dB 

~dt' 



This links the time change of the magnetic field strength to the spatial 
variation of the electric field. The final crucial step is to find an equation for 
the time change of the electric field. From the second equation in (TTOj) one 
may infer that 



V x H = J (for static fields) 

but Maxwell realized that the right hand side of this equation must be 
complemented with the term dD/dt (Maxwell's displacement term), which 
contains the link to the time change of the electric field. This addition is 
needed for maintaining charge conservation (see Eq. (l29l) ). Also, without this 
term no electromagnetic waves would exist in the theory. 

Thus, J C Maxwell was able in 1864 to synthesize the known properties of 
electromagnetism in his now famous equations which give, as far as we know, 
a complete description of the electromagnetic phenomena in the classical 
regime, 



VD 



Q 



V x H = J + 



<9D 

dt 



V x E = - 
V-B = 0. 



dB 

"dt 



Maxwell equations 



(11) 



Here g denotes the charge density and J the current density. We note that 
there is an asymmetry between electric and magnetic fields in the equations 
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in that there appear no magnetic charges (no magnetic monopoles) and no 
magnetic currents. No magnetic charge has been ever discovered, whence 
all magnetic fields are assumed to be generated by moving electric charges 
(electric currents) as described above. 
Using the rule 

V x (V x A) = -V 2 A + V(V ■ A), 

and the relations (J9j) one can show that Maxwell equations give the equa- 
tions, 



_ 2t , 1 d 2 E dJ 1 „ 

VE -^=""»w + ^ (12) 

*H-±^ = -VxJ, (13) 
where we have set (identified with the velocity of light) 

c= 14 
v /eeo7Z/io 

Especially in the case of the empty space (J = 0, q = 0) we obtain the 
wave equation 



1 <9 2 E 

whose plane wave solutions are of the form 



V 2 E-^ = 0, (15) 



E = E cos (k • r ± uot) . (16) 

(Here E is a constant vector.) The magnitude of the wave- vector k 
is 2n/X, where A denotes the wave length, while ui (circular frequency) is 
related to the frequency / by to = 2nf. Inserting (fT6l) into f fT5l ) we infer that 
|k|c = u, which is the same as Xf = c. Eq.(fl6l) describes an oscillating field 
whose frequency is /. The solution can be interpreted as a wave moving in 
the direction of the wave vector =Fk and with the velocity c (light velocity 
in empty space). For empty space we have V • E = which implies, that 
for the plane wave (fTBl ) we must have k • E = 0; that is, the electrical field 
oscillates in a direction orthogonal to the direction of propagation. From 
Maxwell equations we find that the corresponding plane wave solution for 
the magnetic field is then given by 



2 Maxwell equations 



8 



H = H cos (k ■ r ± cut) , (17) 

1 lk 

H = =F k x E = =f- T x E , 

ujHHq r\ k 

where rj = ^/ fi^o/eeo is called the wave impedance (~ 377 Q for vacuum). 
This means that the magnetic field H is orthogonal to both k and the electric 
field En. When both the electric and magnetic fields are orthogonal to the 
wave- vector k the EM- wave is said to be transversal and of the type TEM. 
The direction of the electric field E defines the polarization of the wave. For 
instance, if E = (E x , 0, 0) then the wave is polarized in the x-direction. 




A charge moving with velocity v in an electromagnetic field feels a force 
F given by ((5j). This involves a work per unit time (power P) defined as 
P = F ■ v. Because of the identity v • (v x B) = it follows from (jSJ) that 
P = q~E ■ v. If we have a current density J = gv this result is generalized to 

P= I E-JdV, (18) 
Jv 

defining a power density V = E • J. Using the vector identity 

V ■ (E x H) = H • (V x E) - E • (V x H), 
one can derive from Maxwell equations the following relation, 

E ' J + Ft (l^° H2 + \ ee ° E2 ) + V • (E x H) = 0. (19) 
This can be interpreted as an equation of energy conservation where 

£ = -fi[i H 2 H — ee E 2 
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represents the energy per unit volume associated with the electromagnetic 
field, and 

S = E x H (the "Poynting" vector) (20) 



represents the flow of energy carried away by the electromagnetic radia- 



tion 



39]. That is, given an area element dA then S ■ dA represents the energy 
passing through dA per unit time due to the electromagnetic radiation. 
If we compute the vector S in case of the plane wave (TT6lh (TT71) . we obtain, 

lk k 

S = T-El - cos(k • r ± oot) 2 = ^r/if 2 t cos(k • r ± cot) 2 . (21) 

Tj K K 

Thus we have the important result that the radiated power is proportional 
to the square of the electric and magnetic field amplitudes. Also, for a wave 
which varies as cos(k • r — ut) the power is propagated in the direction of k. 
If we calculate the time average of (ED) over a period T = 2ir/u we obtain 
the factor 



1 1 f T 

— = — I cos(k • r ± ujt) 2 dt. 

2 T J 



Thus, if the average radiation power for a plane wave is 100 mW/m 2 , then 
we obtain the corresponding electrical field amplitude E by setting 

ioo m w/ m 2 = ^|, 



which yields E = \J2 ■ 377 ■ 0.1 V/m » 8.7 V/m (Volt per meter). 

Since electromagnetic waves carry energy, they can also carry "informa- 
tion", which of course make them useful in technology. The motion of charges 
at one place (transmitter) will thus interact with charges at another place 
(the receiver). This interaction is described in terms of the electromagnetic 
(EM) fields. The transmitter generates EM-waves which are intercepted by 
the receiver. 

As seen from the second equation in (fT2l the magnetic field is determined 
by the current J; one can write a solution of the H-equation as 



Here i = t — \r q — r|/c is the retarded time which takes into account that 
it takes time for the field contribution generated at the point r q to reach the 
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point r. This form ([22]) suggests introducing an auxiliary quantity A called 
the vector potential related to E and B by 



B = V x A, (23) 
E=-V#- gr . 

The first equation in (123] ) implies that V • B = holds identically; also it 
implies in conjunction with Maxwell equations that V x (E + ^) = which 
suggests the second equation in (l23l . There is some freedom in choosing A 
and 4>; indeed, if we use the "tilded" versions 

A = A + V£| 

<9£ / "Gauge transformation" (24) 

= -¥ J 

(for some function £) then the fields E and B remain unchanged in (l23l ). 
If we use as the supplementary condition for A that (fixing the gauge to the 
so called Lorenz gauged) 

V.A + I| = 0, (25) 
we get from Maxwell equations 

V 2 A-^^ = -^ J. (26) 
This has a solution of the form (consistent with ([22]) ) 

A(M) = ^/j^A,. (27) 

In principle, if we know the current distribution J (r q , t) in the transmit- 
ting antenna we can calculate the radiated field using fl271) and ([23]) . The 
problem thus reduces to determining the current J (r q , t) which often is a 
very hard problem to solve analytically. However, in many cases simple ap- 
proximations will do quite well. 

Finally we observe two important consequences of the above equations. If 
we differentiate ([25]) with respect to the time and use the second equation in 



2 This was indeed introduced by the Danish physicist Ludvig Lorenz (1829-1891) and 
not by the more famous Dutch physicist Henrik Lorentz (1853-1928) to whom it is often 
attributed. 
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(l23l) together with the first equation in (fill , then we obtain the relativistic 
form of Poisson equation 

o 1 d 2 6 p , . 

*+-*w = ~k- (28) 

This equation determines the potential when the charge distribution g 
is known. A second observation is that the identity V ■ (V x H) = applied 
to the second of the Maxwell equations (1TT1 ) leads to the continuity equation, 



V-J + ^ = 0, (29) 

which expresses the law of the conservation of electric charge. This is of 
importance when e.g. determining the current/charge distributions in anten- 
nas. 



3 Dielectrics and conductors 
3.1 Electric susceptibility 

Common electromagnetic phenomena are due to the interaction of electron 
and protons, the basic elementary particles of ordinary matter. Electromag- 
netic forces hold atoms and molecules together. Since matter is thus made 
up of electrons and protons (and neutrons) we expect matter to affect elec- 
tromagnetic fields and vice versa. Although atoms and molecules may be 
electrically neutral (contain an equal number of electrons and protons), the 
charges may be shifted so that one region is dominantly negative while an- 
other region is dominantly positive. The matter is then said to be polarized. 
The polarization can be understood in terms of electric dipoles. Suppose we 
have a positive charge q at the point r x + 1, and a negative charge — q at the 
point then the potential of the system measured at the point r 2 becomes 

*( \ 1 q 1 q 

0(r 2 ) 



47re |ri + 1 - r 2 | 47re |ri-r 2 | 

When 1 approaches zero such that ql remains a finite vector p (the electric 
dipole moment), the potential becomes 

«r,) = T^ (r = r2 _ ri) . 

47re r d 

The corresponding electric field is given by 
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E(r 2 ) = -V0(r 2 ) = -±- (^Pl - (r = r 2 - ri ). (30) 

The point is that even though the dipole is electrically neutral it generates 
a non-zero electric field which depends the orientation of the dipole. If a 
dipole is placed in an (homogeneous) external electric field E, then a force 
gE acts on the positive end and a force — gE acts on the negative end creating 
a torque N = p x E trying to line up the dipole along the direction of the 
field. This in turn affects the field generated by the dipole. 

The charges associated with dipoles are called bound charges since they 
cannot move freely. Thus, the total charge density can be written as g = 
Abound + £'free- While the displacement field D is defined such that V D = £>f ree , 
one can define the polarization density P such that V ■ P = — Abound- Since 
e V ■ E = g, we obtain 

e V E = V D — V ■ P 

which suggests that 

e E = D-P. (31) 

Experimentally it is found, under not too extreme conditions, that there 
is a linear relation between the external field E and the induced polarization, 
E = e x e P. This finally gives using (13T1) 

D = eo(l + Xe)E = ee E, (32) 

which is equation Q with the relative permittivity given by e = 1 + Xe, 
where Xe is the electric susceptibility. 

3.2 Magnetic susceptibility 

In the magnetic case we do not have magnetic dipoles formed by magnetic 
charges, because, as pointed out earlier, there appears not to exist any mag- 
netic charges in the nature. Instead the magnetic fields are generated entirely 
by electric currents. On the atomic and molecular scales we have electric cur- 
rents due the electrons "circling" around the atoms. Also the "spin" of the 
electrons contribute to magnetism. An external magnetic field may deflect 
the atomic currents and thus change the corresponding field generated by the 
currents in an analogy with the electric polarization. The atomic (bound) 
currents Jbound generate a magnetization M defined by 



V X M = Jbound- 
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a 




Besides the bounded currents Jbound which 
average to zero, we might in conductors have a 
free (macroscopic) current Jf ree related to the 
magnetic field by V x H = Jf ree - Since //q 1 V x 
B = J = Jbound + Jfree we conclude that 

/io X B = H + M. 

For para- and diamagnetic substances the 
magnetization and the magnetic field are, un- 
der "normal" circumstances, linearly related, 
M = x m H, where \ m is the magnetic suscep- 
tibility. Thus, we get B = jj, (l + x m )H = /x /xH, with the magnetic perme- 
ability given by \i = 1 + Xm, which is the second equation in (J9j) . Whereas 
the electric polarization was analyzed in terms of electric dipoles, the mag- 
netization may be analyzed in terms of small current loops. If one consider 
such a small current loop in a magnetic field characterized by B, then the 
total force acting on it is zero whereas the torque becomes (applying ((5|)), 

N = j) Ir x (dr x B) = m x B, 
with the magnetic moment m defined by m = Ja, where 



1 



r x dr 



is the area enclosed by the loop and I its current. The magnetization M 
corresponds to the density of magnetic moments. 



3.3 Ohm's law 

A "free" charge q in an electric field E feels a force gE which causes it to 
move. Thus, electrons in the conduction band in metals (conductors) can 
form a current when a potential difference is applied over a piece of a metal. 
The electrons though meet resistance caused e.g. by the thermal motion of 
the atoms. This is manifested in the well known "law" of Ohm according to 
which one needs a potential difference U = RI in order to drive a current I 
through a conductor with the resistance R. (It is conventional to use U for 
the potential in the theory of circuits.) In terms of the current density J and 
the electrical field E driving the current, the law of Ohm can be written as 



J = crE 



(33) 
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where the conductivity a is inversely related to the resistance. More 
precisely, for a conductor of length L and cross section A we have for the 
resistance 

L L 

R = — r = Pt- 

where p—\ja defines the resistivity (typically of the order of 10~ 7 Qui 
for metals). When treating electromagnetic waves in conductors we thus 
have to use the relation (1331) in Maxwell equation (fTTi ). 



3.4 Shielding 

3.4.1 Skin effect 

It is well known that metallic enclosures ("Faraday cages") protect against 
external electromagnetic fields. This shielding is caused by the fact that the 
charge carriers generate an opposing field. Suppose we have an incident plane 
wave along the z-direction on a metallic surface with the normal direction in 
the — ^-direction. The equation f fl"2l ) becomes^] 

„,„ 1 d 2 E <9E . , 

VE -«-*a =a (34) 

Inserting a solution of the form (in the metallic medium z > 0; here 
we consider only the transmitted component, not the incident and reflected 
components in z < 0) 

E x = E 0x e^ kz ~^ 
we obtain for k the equation 

k 2 = — + ianoficu. (35) 
c 2 

The imaginary part leads to exponentially decaying factor in exp(ikz). 
For instance, if the second term in (1351) dominates then we have 



/ 1 + t / 

k « —^-y/afiofiu, 



3 The p-term vanishes in this special case. Combining the continuity equation f29|) with 
J = crE we obtain for an harmonic plane field parallel with surface, E = (E x ,0, 0), 

q = -«-V • E, 

to 

which is since E x depends only on z as the plane field travels in the z-direction. 
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which leads to 

e lkz « exp (i—j^y/a ^q^jlujz^ ■ exp ^— — y^y/a^Qfiuz 

The second decay factor shows that we have characteristic penetration 
depth of 



S = \ . (Penetration depth, for high a) (36) 

As an example, for copper we have a = 6 ■ 10 7 (fhn) -1 , /x ~ 1, which give 
at uj = 2 7i"/ = 2 7r 2.4 GHz a penetration depth of 5 = 1.3 • 10~ 6 m; that is, 
about 1 /im. From this it follows that a 1 mm Cu-sheet will practically stop 
the field completely. For an aluminum (5 = 3.8 • 10~ 6 m) foil of thickness 0.01 
mm we would get a suppression factor around exp(— 0.01 mm/5) m 2.5 • 10~ 3 , 
or -52 dB. 

The shielding (absorption of radiation) can also be interpreted in terms 
of a complex permittivity. From the second of Maxwell equations (TTT1) we see 
that the current density J and the field D occur in form of the combination 

Assuming harmonic fields depending on time as exp(— iujt) the time deriva- 
tive above can be replaced by the factor — iu, and if we replace J by crE, the 
above expression becomes 

<9D 



^ (JL> /a \ „ 

crE + — - = -iu i— + ee E. 

at \ uj J 



This means that for conductors the effect of the current on the fields can 
be taken into account by using a complex relative permittivity given by 

e r = e + i = e + ie" . (37) 

iue 

It is conventional to denote the real part by e' and the imaginary part 
by e". In many texts they write the complex permittivity as e' — ie" , which 
follows from assuming a time dependence exp(iut) instead of exp(— icut); 
thus, it is purely a matter of convention. The effect of the conductor is 
also to make the wave impedance r] defined in (TT7T) imaginary. Indeed, using 
Maxwell equations for the plane waves we find that the E- and H-amplitudes 
are related by 
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where k is given by (I351) . The above relation reduces to the one in jlj 
when a = 0. A material is called a good conductor if the dielectric imaginary 
part dominates, e" ^> e', which translates into 

^> 1. (Good conductor criterion) 

ee uj 

As an example, for uj = 2 nf = 2 n 2.4 GHz we get 

e oo « 0.133 (Vim)- 1 (f = 2.4 GHz). 

This can be compared with the conductivity of copper, a = 5.8 -10 7 
(f2m) _1 , which thus, as all metals, qualifies as a "good conductor" by a safe 
margin. The human body has a conductivity around 0.2 (f2m) _1 and is 
thus a poor conductor at this frequency. Sea water is a borderline 
this frequency having a conductivity around 4 (fim)" 1 and e ~ 80. The 
penetration depth is an important characteristic for food that is heated in 
microwave ovens; 5 must be of the order of centimeters for the radiation to 
heat the food thoroughly. 

Because of the small penetration depth for good conductors it is typically 
assumed that the electric field is zero in the main part of the conductor, and 
that it is, like the current, confined to a thin layer of thickness about S near 
the surface. This has important consequences. Suppose the surface lies along 
the xj/-plane, and that the electric field is tangential in the x-direction. Then 
we have from the third equation in ffTTjl . 

dBy 



I AS- 1 = I Ad 



Of 



If take the difference AE X to be over the interface, and let |Az| — > 0, then 
we obtain that 



Ef) = E«, (39) 

that is, the tangential component E t of the electric field changes contin- 
uously across an interface between two mediums (I) and (II). Thus, if the 
electric field is zero inside the conductor (e| 7 ^ = 0) it must also be zero at the 
outside surface (e[ // ' ) = 0). It follows that when an oscillating electric field 
E l impinges on a conducting surface, it generates a surface current J causing 
magnetic field H which in turn, according to Maxwell equations, causes a 
reflected field E r , such that the tangent component of the total field is zero 
at the surface (this is a simplification valid only on a scale large compared 
with the skin depth) 
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+ W t = 0. (At the surface of a conductor.) 



(40) 



This is used as a boundary condition when treating radiation in cavities 
and antenna radiation. Thus, the current J caused by an impinging field in 
an antenna will generate a field outside the antenna from which one may, for 
example, calculate the gapfield and corresponding potential difference in case 
of a dipole antenna (see sec. 15 . 2 . 1 B . Another consequence of the penetration 
layer is that electromagnetic waves will loose energy due to ohmic losses; the 
waves induce surface currents which encounter a resistance given by 



Using previous data on copper we get for its surface resistance R s ~ 1/78 
tt (at 2.4 GHz). 

3.4.2 Leakage through slots 

If there are holes in the shield then there will be no counteracting current 
at that place and radiation can leak through. This effect can be used to an 
advantage when constructing slot antennas, but for shielding purposes the 
leakage is of course a nuisance. In order to get an estimate of the leakage 
through a slot one may consider a rectangular hole in a conducting sheet. 
We suppose the sheet is in the xy-plane with the normal along the — z-axis 
and has thickness d. Further we take the hole to have the corners (0,0), (0,6), 
(a,0), (a, b) in the xy-plane. 

We will treat the hole as a wave guide on which impinges a planar EM- 
wave along the z-axis. Intuitively it seems clear that waves with a wavelength 
A ^> a, b will have difficulty in passing through the hole; that is, the hole acts 
as a high-pass filter damping waves with wavelengths exceeding the dimension 
of the hole, but letting smaller wavelengths through (the high frequency 
part). We will consider a TE-wave passing the wave guide; thus, the electric 
field is transversal while the magnetic field may also have a longitudinal 
component along the z-axiifl We will therefore assume that 

4 The TEM case leads to a trivial solution of zero fields inside an empty wave guide. 
Indeed, in this special case one obtains from Maxwell equation that (T means here that 
the operators are restricted to the transversal xy-plane) 




(Surface resistance.) 



(41) 



V T x E = 



from which one may posit that there is a function 4> such that 



E = -V T (/>- 
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E(x, y, z, t) = (E x (x, y), E y (x, y), 0) • e^ z ~^\ 

H(x, y, z, t) = (H x (x, y), H y (x, y), H z (x, y)) ■ e^ z ~^ . 

The dependence on z is thus factored out as exp(iflz) since we are inter- 
ested in wave solutions progressing in the ^-direction. If we insert the above 
ansatz into the wave-equation ( fT2i ) for the E-field we obtain (J = and g = 
in empty space) 



V%E X 



T ±-i x — ((3 — k )E X , 
with 

, d 2 d 2 



(Helmholz equation) 



(42) 



T 



dx 2 dy 2 



and k 



with a similar equation for the y- 
component. Here V? refers to the Lapla- 
cian operator restricted to the transversal 
plane. Since the tangential components 
of the electric fields vanish at the surface 
of the walls, they will be of the form 



E x (x,y) 
E y (x,y) 



g[x) sm j , 
h(y) sin (^H^ 



Using the Maxwell equation d x E x + 
d y Ey = we can determine the functions 
h and g, obtaining finally (m, n ^ 0), 



10 

c 



2tt 



\ b 

X 






z 








a 



E x (x,y) 
E y {x,y) 



Combining this with the equation Vt x E = gives the equation Vf^ = 0. This is the 
Laplace equation in two dimensions, and for a simple region (such as the cross-section of 
the wave guide) this has the trivial solution </> = constant. Indeed, -©tangent = implies 
that 4> is constant along the boundary of the cross-section, and therefore that Vy</> = 
has the trivial solution <fi = constant in the cross-section. In the TE case, where we allow 
for a magnetic longitudinal component, Vt x E = is replaced by Vt x E = iu>fiQ^,H z , 
and non-trivial solutions become possible. 



E mn -cos[-— )**{— J, (43) 

_, b fm%x\ /n7iy\ 
-E mn — sm cos — — . 



■"mn 

n \ a 
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The case m = corresponds to a solution of the form 



E x (x,y) 
E y (x,y) 

The general solutions may be constructed as superpositions of these (m, n)- 
mode solutions. Inserting an (m, n)-mode solution into (|42l) we obtain the 
relation 




From this we see that if a, b < A then (3 must be necessarily imaginary 
leading to an exponential decay factor exp(i/3z). Suppose we have a narrow 
slot a = 1 mm in a d = 2 mm thick conducting plate, then we may estimate 
the radiation through the slot to be damped by factor of the order 



. fnny\ 
anSm {~b~) ' 
, . (rmx\ 
bn*n{—). 




exp<-\ (-) -(-r) • d } w exp(-3141 • 0.002) w 0.0019, 



which corresponds to a -54 dB damping (power) at 2.4 GHz (A = 0.125 
m). For a < A the above formula for damping can be approximated by (in 
dB) 



damping [dB] = 20 • log (e~^) ~ -27.3 



d 



3.5 Reflexion and refraction 

Reflexion and refraction of waves is a fa- 
miliar phenomenon from our daily expe- (I) (II) 
riences with light, sound and water. Re- 
flexion is a basic effect when a wave hits V/////////// 

an inhomogeneity in the medium, typi- V////////// 

cally the interface between two different 

mediums such as air and water. We will 

discuss this fundamental feature in terms 

of a very simple model. Consider a "wave" 

traveling along the x-axis in a medium (I, 
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x < 0) whose propagation velocity is U\. At x = starts another medium 
(II, x > 0) with a different propagation velocity u 2 . (In case of EM-fields 
u = c/a/cco where c is the light velocity in vacuum.) We use for a propa- 
gating field satisfying the following equations: 



d 2 (h 1 d 2 (h 

d-^w = ' (I)(I<0) (45) 

** l8 "* 0. (II)(x>0) 



dx 2 u\ dt 2 



We write the basic harmonic solutions for regions (I) and (II) as: 



(j) 2 (x,t) = Te i{k2X ~ wt) . 

Here the wavenumbers ki are given by fcj = 
uj/ui. The interpretation of the solution (1361) is 
the i?-term represent the reflected part prop- 
agating in the — x-direction, and the T-term 
the transmitted part. The coefficients R,T 
can be determined from the requirement that 
4> and d(f)/dx be continuous at the boundary x 
= 0; that is, 0i(O-) = 2 (O+), d</>i(0-)/&r = 
•902(0+) /dx. This yields the equations 



1 + R = T, 
h(l-R) = k 2 T, 



(47) 
(48) 



(I) 
(II) 



(46) 




from which obtain R and T, 



i2 = ^ = ^l, (49) 

2fci 2« 2 , N 

T = V- = — • 50 

fcl + K 2 Ml + U 2 

The above model may e.g. be used to describe the effect of connecting two 
cables with different impedances; the discontinuity at the connection gives 
rise to reflexions. We may also note that the sign of the reflexion coefficient 
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R depends on whether the wave travels faster or slower in region II than in 
region I. 

Next we will consider an EM plane wave in open space (z > 0) impinging 
on a dielectric surface in the xy-plane at z = 0. The total electric field on the 
side I (z > 0) will consist of the incoming and the reflected part (R), while 
on the side II (z < 0) we will have the transmitted (refracted) part (T), 



Ee i(k.r- W t) +E Ji e i(k«.r- W t) j (i) (5i) 

E T e i(k r .r-^ (H) 

Here the magnitudes of the wavevectors are given by 



ki = |k| = |k | = Uy/fitfioeteo, 
k 2 = \k T \ = LOy/fi 2 fi e 2 eo. 

As demonstrated in connection with (j39j) the tangential component of 
the electric field does not change across the interface. Hence, choosing the 
coordinate system so that the tangent component is along the y-axis we 
obtain at the I-II interface, z = 0, 

This equality is only possible if (at z = 0) 

k ■ r = k R ■ r = k T • r 

from which one deduces (set r = y ) that 
the angle of reflexion 9 r is equal to the inci- 
dent angle 6 (angles are here measured as those 
made by the directions of propagation with the 
interface), while the angle of refraction 6 t on 
the other hand is related by "Snellius' law" 

rii cos 6 = n 2 cos 9 t (Snellius) (53) 

where n« are the indexes of refraction of the 
mediums given by 7ij = ■Je^i. From (T52T) we 
also obtain that 




E y + E* = El 



(54) 



3 Dielectrics and conductors 



22 



If we consider the .E-field to be polarized in the yz-plane (vertical, V- 
polarization case; H-field will be along the y-axis) then the components of 
the fields can be expressed in terms of the amplitudes, 



Ey 


= E sin 0, 


Hy 


= o, 


E z 


= E cos 9, 


H x 


— —E, 






E« 


= -E R sin9, 


»»" 


= 0, 




= E R cos e, 




= -E R 




= E T sin9 t , 




= 0, 




= E T cos9 t , 




V2 



In order to determine the amplitudes E R , E T , we need one further equa- 
tion besides (i52l ). This can be found by applying f TTUT ) to a very thin "pill-box" 
which contains the I-II interface, then one obtains that 

£n-dS=(D I z -D I z I )S = Q, 

where S is the top (bottom) surface area of the pill-box, and Q the surface 
charge contained by it. Going to the infinitesimal thin pill-box limit we obtain 
the general result on the normal component of the displacement vector, 

= + q„ (56) 

where g s is the surface charge density (Q/S). In our particular case we 
can assume that there are no extra surface charges (g s = 0), whence, using 
D = ee E, 

e 1 (E z + E R ) = e 2 Ej. (57) 

Another alternative is to use the boundary condition that the tangent 
component of the H-field is continuous across the interface, Hf = H^, which 
can be derived in a similar way as in the case of the E-field. Note that 
the wave impedances rji in (1551 ) are given by ( 1381 ) which cover the case of 
conductive media too. 

Combining (l5Tj) . (jSD and (l53l) . we can obtain after some algebra, 
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E R e r smO-y/e r ~ (cosfl) 2 

& e r smv + We r — (cosvy 
_ _ E T _ 2 sin 6> ^/e; 

E e r sin 9 + a/ e r — (cos 6>) 2 ' 
(V-polarization case.) 

where we have used the notation e r = e 2 /ei. Similar considerations can 
be applied in the horizontal (H) polarization case, with the end result, 



E R _ sin 9 - y/e r - ( 



cos I 



\2 



& sm 9 + A/e r — (cos 0) 2 
E r 2sin0 



sin + \/e r — (cos9) 2 
(H-polarization case.) 

The equations (1581 ). ( |59l) . are known as Fresnel [FRA-nel] equations. 
These equations also apply when interfacing a conducting material by re- 
placing e by a complex number as explained in connection with (1371) . Thus, 
if the medium II is a perfect conductor this corresponds to letting \e r \ — > oo, 
and we have then p v — 1 for V-polarization and ph = — 1 for H-polarization. 
For dielectrics there is a special angle, the Brewster angle 9b, 

sm9 B = -=L= (60) 
Vl + e r 

at which the reflected vertical component goes to zero, p v = 0. Thus, if 
the incoming field is vertically polarized, none of it will be reflected at the 
Brewster angle (for a planar interface) . This means that if the transmitter 
(using vertical antenna) and the receiver are placed such that the Brewster 
angle condition is satisfied then the reflecting component of the radiation is 
eliminated from the transmission, and only the direct field is received. This 
configuration may be used to measure how the inclination of the receiver 
antenna affects the reception; that is, to measure the function G{9) for vary- 
ing 9. The above analysis can be generalized to the case where we have 
two mediums I, III, with a second medium II of thickness d sliced between 
them. One example might be air (I), and ice sheet (II) with water (III) below 
(which we have investigated experimentally). The reflected field in I is thus 
reflected both from the interface I-II and from the interface II-III. We may 
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treat the problem with the methods used above, pasting together plane wave 
solutions in the regions at the interfaces. We may also use the methods of 
geometrical optics and sum the contributions from all the additional reflex- 
ions from the intermediary layer interface II-III. We denote by (and similarily 

for transmission coefficient) p (j^, O^j the reflexion coefficient (subindexes will 

indicate the polarization states) for an EM-wave in a medium I impinging 
on a surface of a medium II at the angle 6. 

The total reflexion will be a sum of the pri- 
mary reflection at point A (see figure), the next 
contribution comes from the transmitted part 
which reflects from point B and then exits the 
surface at point C, and so on. It is important 
to note that the parts that bounce through the 
intermediary layer II pick up additional phase 
differences due to the factor 

exp(ik n ■ r). 

The phase contribution due to a a given optical 
path is ks where s is the length of the path. 
Thus, after some trigonometrical exercises, the 
phase difference between the paths ABC and AD will turn out to be, 

n 2 27i2d rii2-K , „ „ 2n 2 2ixd ( 1 cos#?\ 

A = -±— - !— 2dcot0 t cos0 = —\ — -T- - — -f = 2kn 2 dsm6 t , 

Asm 6^ A A \smu t smv t J 

where we have used the law of Snellius, and k = 2n/\ for the wavevector 
magnitude in vacuum, and rii for the indexes of refraction. Summing all the 
reflexion contributions we get, 




Aotai = P (|, 6j + e iA r (|, 6j p (|, t ) t (|, d)j + (61) 

e M^MH 4*0 + ■ - 
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Here we have used the geometric summa- 
tion rule 1 + x + x 2 + ■ ■ ■ = 1/(1 — x). By a 
similar calculation we obtain for the transmis- 
sion coefficient r for the radiation that enters 
into the medium III, 




r = 



e^/V(^)r(gA) 



) 



(62) 



B 



As an example we can calculate the trans- 



mission coefficient for radiation impinging nor- 
mally on a brick wall of thickness d = 10 cm 
assuming e\ = £3 = 1 and 62 = 4, leading to (setting e r = ^/ti) 



r = ^— ? = -0.74 + z0.43, 

(1 + + e 2k ^ (^r r - iy 



4 nr e ik^e^d 



(63) 



whose absolute value is 0.86 (@ 2.4 GHz) corresponding to a reduction 
of power by the factor 0.86 2 = 0.74 (-1.3 dB). Since e r is assumed to be real 
there is no absorption in the wall. If we use e" = 0.07 for the brick wall then 
we have for aim wall |r| = 0.36 thus showing already significant absorption 
(-8.9 dB). In reality there would be a further loss of power due to scattering 
caused by the inhomogeneities in the wall. Brick walls are seldom 1 m thick, 
instead the radiation may have to pass several brick walls which are say 10 
cm thick. Then a quick estimate would be that the power decreases with a 
factor about 0.74 per wall. For a more detailed treatment one can extend 
the above methods to an arbitrary numbers of dielectric layers |40j|. One can 
also apply the transmission theory for cables using for impedance the wave 
impedances Z = a/ yU/i /eeo (see Appendix ID]) . 

3.6 Image charges 

We consider a charge q above a plane perfect conductor which we take to 
be in the xy-plane at z = 0. The plane will now affect the electric field. As 
explained earlier the electric field must have a zero tangential component on 
the surface of the conductor. An equivalent formulation is that the potential 
<f is constant along the surface. Thus, given this boundary condition, one 
has to solve the Laplace equation V 2 <p = which is valid for z > 0, except at 
the place of the charge which we may suppose is at the point r = (0, 0, h). 
One can convince oneself that the solution must be 
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<p(r) 



Q 



r 



Q 

47ree r 



r + r 



(64) 



r |* 47ree |r + r | 2 

It satisfies the Laplace equation for z > except at the point r , and it 
vanishes for z = 0; that is, z = is an equipotential surface. Furthermore, 
if we integrate E = — V<£> over a small sphere containing q we obtain q/ee 
proving that it is indeed the potential of the charge q. The solution (1641) 
means that the effect of the conducting plane is the same as if we had an 
additional extra charge of the opposite sign at the place of its mirror image, 
r — 2(n ■ r ) = — r (n is the normal of the surface), in an empty space. 
One consequence of the mirror effect is that the charge is attracted toward 
the conducting plane by the apparent opposite image charge. Physically 
the effect of the conducting plane is that the charge q polarizes the free 
charges in the plane by attracting them if of opposite sign, and repelling them 
otherwise. In fact, the induced surface charge at z = can be calculated from 
g s = —ee dip/dz by inserting the solution (1641) . If we integrate g s over the 
surface we get in fact for the total induced charge the result —q. 

This mirroring method can be generalized to other surface that can be 
construed as equipotential surfaces for some distribution of charges. Consider 
the case where we have two conducting planes meeting along the z-axis. We 
may take the conducting planes to be the xz-plane and the yz-plane (see 
part (b) in the figure). The potential in the open space region is obtained 
adding three image charges as shown in the figure. 



(I) i q" 




r - r _ f 



-Q 



////////// 



(b) 



Indeed, one sees that this arrangement makes the total potential zero 
along the planes. This example has application in the case we use 90°-degree 
bent sheet as an antenna reflector. A somewhat more involved case is that 
of placing a charge between to parallel conducting planes which requires an 
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infinite number of image charges (as one can "see" from the analogous case 
of placing a candle between two parallel mirrors). 

The imaging principle can also be applied to the case of a dielectric instead 
of a conducting plane. Thus consider the case where we have a homogeneous 
dielectric I in the region z > 0, and a different dielectric II in the region 
z < 0. We place a charge q at the point r = (0, 0, h) in I, and the problem is 
to determine the resulting potential in I and II. We make the ansatz that the 
potential in I is the sum of the potential generated by q and an imaginary 
charge q" in II, and that the potential in II is generated by a charge q" at r 
possible different from q due to "screening". 

^ ^ ~ 47re ieo |r-r | + 4vr ei e |r + r |' ^ ^ ^ 

/7 (r) = /" 1 (In II.) 

47re 2 e |r — r | 

We can determine the unknown charges q' and q" from the boundary 
conditions at z = where we have the continuity of the tangent electric field, 
E( = Ej J , and the normal electric displacement, = (since no free 
surface charges are expected for the dielectrics). Expressing these conditions 
in terms of the potential (l65l) we obtain the equations 

*+^=^ (66) 
ei ei e 2 

-q + q' = -q". 
These equations have the solution 



q = q — : — , (67) 
ei + £2 



2e 2 

q = q- 



ei + ^2 

The conducting plane corresponds to the limit e% — ► oo and we see that 
in this case we indeed recover the solution q' = —q. 



4 Interference and diffraction 
4.1 Geometric optics 

As is well known, light, which is EM-wave of very short wavelengths (around 
0.1-1 /im), can in many problems be treated as consisting of "rays". This is 
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the method of geometric optics. The plane waves discussed above do not exist 
in reality as they would be of infinite extent. But locally EM- waves may often 
quite well be approximated by plane waves. The geometric optical methods 
give good approximations far away from the source (transmitter) and when 
we consider spatial dimension large in comparison with the wavelength A. 

If we write if) to represent some component 
of the EM-field, then one may set 



if{r,t) = A-e iS{r ' t \ (68) 

where A (amplitude) and S (phase) are real 
functions. In dielectric medium if) satisfies the 
wave equation 



v 2 v> 



n 2 d 2 if) 







(69) 




where n = *JejI is the index of refraction. 
A plane wave corresponds to the case where A = constant and S(r, t) 
k • r — cut. For a fixed time t the sets of points r satisfying 



S(r, t) = constant 

form surfaces (wavefronts) of constant phase. If we draw the lines that 
are everywhere orthogonal to the wavefronts we obtain the rays. The concept 
of rays is useful only when the characteristic dimensions of the regions con- 
sidered are large in comparison with the wavelength A. Then one make the 
assumption that the amplitude A changes only a little over regions of the di- 
mension A. Mathematically this condition may be expressed as |VA|A <C \A\. 
If the phase locally is close to that of the plane wave we may assume that 
| V*S| A ~ 27T. Using these assumptions one obtain from f j68~l) the approximate 
equation 



n k . (Eikonal equation.) 



(70) 



The last equality follows if we assume an harmonic wave for which S(r, t) 
is of the form Sq(t)— ut. In (1701 k = 2n/X is the magnitude of the wave vector 
in vacuum. The eikonal equation is equal to Fermat's principle^ according 
to which the rays are paths T which minimizes the traveling time defined by 



5 This follows from the fact the eikonal equation is analogous to the so called Hamilton- 
Jacobi equation of classical mechanics for a particle moving in a potential proportional to 
— n(r) 2 . This in turn is related to the least action principle (Maupertuis) which finally 
leads to the Fermat's principle. 
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f nds 

tr= • 

Jr c 

From this it follows that rays are straight 
lines in homogeneous regions (n independent of 
position). If we are far away from an antenna 
(r ^> A) then we may think of the EM radia- 
tion as arriving from the antenna along rays. If 
the antenna is close (in relation to the distance 
to the observation point) to the ground, build- 
ings etc, then, besides the contribution coming 
along the straight line between the antenna and 
the observation point, we may also have rays 
reflected from the ground, buildings etc, arriv- 
ing at the observation point. Thus, the field at 
a point r may be a sum of contributions due to 

^(r) = A ie iSl(r) + A 2 e iS2{r) + ■■■ . (72) 

For a flat ground this sum reduces to just two parts ("two-ray model"): the 
direct contribution and the reflection from the ground. In case of reflexions 
the terms in ([T2l will be of the simple form Aexp(iks) where s is the total 
path length. 

Because the different paths may have different lengths Si the sum ( 1721) 
may lead to either constructive or destructive interference. We have maxi- 
mum destructive interference if the path difference is A/2 + nX (phase dif- 
ference 180°) and maximum constructive interference if the path difference 
is n\ (phase difference 0°) in terms of the wave length A. As a simple ex- 
ample we consider the interference between the direct and reflected ray from 
a transmitter (T) to a receiver (R), both at the height h above a perfectly 
conducting ground. When the distance r between R and T is considerably 
larger than h we may assume that A\ ~ A 2 and the magnitude of the in- 
terference becomes proportional to (this corresponds to the case of vertically 
polarized EM waves) 

gifcr _|_ gifc2^/ h' 2 +r 2 /A 

This leads to a characteristic interference pattern with changing distance 
r with a separation Ar between the positions of maximum amplitude ap- 
proximately given by 



(71) 



T r R 

FT — — *~ — — ~7\ 

\ / 
h Ns * 

//////////?/////////// 



many paths, 



1 + e 



ikh 2 /r 



4 < cos 



7T 



rr 



Ar 



(73) 



4 Interference and diffraction 



30 



For r > 2h 2 / A the cosines term approaches 
1. With reference to the Fresnel equations (1551 ). 
(1591) , we make the observation that, for a re- 
flexion from a perfectly conducting surface, the 
reflexion coefficient ph is negative for the hor- 
izontal polarization case. This is also true for 
the vertical polarization reflexion coefficient p v 
in case of a dielectric ground when the distance 
r (see Eq. (I158f) ) is large enough to make 9 
smaller than the Brewster angle. In these cases 
the "+"-sign in (l73l) must be replaced by a "— "- 
sign, and the cosines term becomes instead a 
sinus-term, 




This has the interesting property of approaching 




as r — >• oo, meaning that the interference reduces the power (which is propor- 
tional to the square of the amplitude) with an additional r~ 2 -factor. Since 
the amplitude A falls of as r -1 the power will fall off as A 2 r~ 2 oc r -4 in this 
case. Thus the long-range behaviour is vastly different depending on whether 
we have a summation or a subtraction in (|7"3l) . 

One way to understand the propagation of the wavefronts is to imagine 
that every point on the wavefront is the source of an expanding spherical 
wavefront (a "Huygens wavelet"), which together with the other spherical 
wavefronts form the new wavefront. This principle was advanced by Huy- 
gens (1678), and it gives a nice explanation of why reflection angle is equal 
the incidence angle, and for the law of Snellius (noting that the velocity of 
propagation of the wavefronts is c/n where n is the index of refraction). It 
also provides a picture of the diffraction of EM waves e.g. through a hole in 
a screen. 




4.2 Fraunhofer and Fresnel diffractions 

Specifically consider a rectangular hole extending from —a/2 to a/2 in the 
x-direction and from y = —6/2 to y = 6/2 in the ^/-direction of a an opaque 
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screen. We suppose that a plane wave ip propagates in the ^-direction and 
impinges on the hole at z = 0. The field at the observation point r = 
(zo, yo, Xq), will be obtained by summing the phase factors exp(iks) over rays 
from the surface of the hole to the observation point 

/a/2 r-b/2 
I e iks{x ' y) dxdy. (74) 
-a/2 J-b/2 

If we take r to be the distance from the center of the hole to the obser- 
vation r point then the path length s(x,y) can be evaluated as 

L , f x-xq\ 2 f y - yp \ 2 (x - xq) 2 (y - y f 

S = r f + {—) + {—) ~ r + —2^ + —2^- (?5) 

Here we have assumed that r ^> \a\, \b\. If x ^> a and y 3> b (Fraunhofer 
case) one can retain only the linear term in (x — Xo) 2 = x 2 — 2xxo + x^ (and 
similarly for y) which simplifies the integrals (1741 ) to the form 

r f 2 e-^-^dxdy = sin ( !^) 2L sin ( ^) . 
J -a/2 J-b/2 kx \ 2r J ky \ 2r J 

This means that the intensity at the point of observation will be propor- 
tional to the factor 

/ sin(fcax /2r) \ 2 / sm(kay /2r) \ 2 
V kx /r J V ky /r J 

Closer to the hole we may no longer have 
Xq ^ a and y ^> b. In this case (Fresnel 
diffraction) we have to retain the full quadratic 
expression, which leads to integrals of the form 

/•a/2 

/ e lk< ^ x ~ X(> > l 2r dx 

J -a/2 

This integral cannot be solved analytically 
but is easily evaluated numerically using com- 
puters. Traditionally the integral has been an- 
alyzed in terms of the Fresnel integrals 
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C(a) = / cos(7rs72)ds, 
Jo 

S(<r) = / sin(vrs 2 /2)rfs, 



'0 

and the so called Cornu spiral defined by x = C(er), y = S(u). 

We may consider the example of diffraction by an edge at y = in the 
xy-plane with the open space for y > 0, and the screen in the half y < 
0. Suppose the observation point is at r = (0, d, D p ); that is, D p is the 
distance from the xy-p\&ne to the observation point. We have the geometric 
shadow for d < 0, and the illuminated region for d > 0. Repeating the same 
arguments as above the calculation of the field at r leads to an expression 
involving the integral 



o 



oo /on r°° 

e*Wl w >dy = \ / e^d-q, (77) 



where the diffraction parameter w is defined by 



We consider the shadow region d < 0. Then it is possible to estimate the 



integral ( 1771 ) by partial integration [26|, p. 152] - doing it twice yields 



e^d V = e™ 2 -— - + — \-- v -^d V (79) 

For large \w\ the first term will dominate and the intensity of the field \ip\ 2 
will therefore be proportional to l/4w 2 . Comparing with the corresponding 
calculation for d — >• oo (far into the illuminated regionj^, it follows that the 
intensity P in the shadow region varies as 



(80) 



Attw 2 

where Pq is the intensity in the illuminated region. We may apply this 
to an example discussed in joi, p. 132] and treated by other methods there. 



6 For this it is useful to know that 

exp(iri 2 )dri = y/n/2 ■ (1 + i). 
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Thus consider a transmitter at the height of 12 meter and a receiver 37.3 
meters away at the height 2 m. Between them is a house with 12 m to the 
notch, which in turn is 20 m from the transmitter. For a rough estimate of 
the reduction of signal strength due to the building, which is assumed to act 
as a screen, we can use (1801 and ( 1751 ) with D p ~ 17 m and d = -10 m. In the 
case of / = 900 MHz we find the reduction to be about -28 dB, which can be 
compared with the value -25.8 dB given by Bertoni. From (1751) we observe 
that the boundary of the shadow region w = constant is given by the lower 
part (y < 0) of the parabola 



2 v 2 - (si) 



k 

z = — 
w 

Above we assumed that the waves from the transmitter are plane waves 
(parallel rays) , but if we take into account that it is a finite distance D q from 
the screen we need only to modify the expression for diffraction parameter 
w ( 1751 ) according to 



w = d 



kD q 



2D P (D P + D q 



The Fresnel and Fraunhofer cases above be roughly distinguished by that 
the Fresnel approach deals with the situation where we have diverging rays 
while the Fraunhofer approximation treats the case with parallel rays. Since 
the Fraunhofer approximation uses linear terms in the argument of the exp- 
function it lends itself readily to Fourier methods, and this has lead to the 
development of Fourier optics. 



4.3 Fresnel zones 

Consider a transmitter (T) and receiver (R) whose line of sight (LOS) dis- 
tance is D = d\ + g?2- Let U be a point on the line of sight with TU = c?i, 
UR = d 2 , and let UV be perpendicular line to the line of sight with length s. 
Denoting Si = TV and s 2 = VR the difference between lengths of the paths 
TR and TVR becomes 



F = si + s 2 - di - d 2 = \Jd\ + d 2 + yj d\ + d 2 - d x - d 2 . (83) 

The n:th Fresnel zone is the region bounded by F = n\/2 which forms an 
ellipsoid with the transmitter and receiver at the focuses. The first Fresnel 
zone Fi thus corresponds to maximum phase shift of 180° due to a reflection. 
As a rule an undisturbed path between the transmitter and receiver requires 
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that there be no obstacles in the Fresnel zone F x . If D ^> d then one can, 
for d\ = d 2 = D/2, approximate (|83l) by 

f -Hf 

D ' 

from which we obtain the width of the n:th Fresnel zone, 

rn=d= y— = y—- 

Inserting this into (1751) and (1501) we find that the reduction of the intensity 
of the field due to a shadow of depth r n would be 1/itn. Roughly, if the 
Fresnel zone F n is free then the transmission power is disturbed at most by 
a factor of (1 — l/im). In decibels 1/tt corresponds to 10 log(l/7r) ~ -5 dB. 
The parameter 

d_ 

r n 

where d is shortest distance from LOS 
to obstructing objects is called the Fres- 
nel zone clearance and is used in investi- 
gating the path clearance. While objects 
protruding into the Fresnel zones cause 
scatterings, one draws a series of Fresnel 
zones also in order to see whether there 
are some surfaces close to the tangent 
of the Fresnel zones; such surfaces may 
cause reflexions and thus interference ef- 
fects. 




4.4 Kirchoff equation 

Above we presented the Fresnel-Huygens principle in an intuitive fashion. 
It can also be justified starting from the wave-equation (1691 ) which in the 
harmonic case ip(r,t) = ip(r) exp(— iut) becomes, 

V 2 ^(r) + fc 2 ^(r) = 0. (k — uj/c.) (84) 

Suppose ip(r) is known on some surface St, the idea is to try to express 
the value of ip(r) at an observation point r in terms of its values at the 
surface S%. First we note that the function (a Green function) 
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G(r) = -^— e ikr (85) 



satisfies 



V 2 G(r) + k 2 G(r) = 6(r). (86) 
From that follows (using Green's theorem) 

I (^(r)VG(r - r ) - G(r - r )V^(r)) ■ dS = (87) 

JdV 

[ (^(r)V 2 G(r - r ) - G(r - r ) V 2 ^(r))dV = 
Jv 

(^(r)(V 2 + k 2 )G(r - r ) - G(t - r )(V 2 + k 2 )^{v))d\/ = 
ip(r)8(r - r )dV = ip(r ). 



We can take V to be a volume enclosed by a screen Si (and some hole in 
it) and an infinitely far away surface S2 where the field disappears. It may 
assumed that the field is zero at the screen and is approximately a plane wave 
at the hole area A, where we thus have Vip = ik in -0 . Hence lf87]l becomes 
(Kirchoff equation, 1882) 

V>(r ) = / ^(r)(VG(r - r ) - tk in G(r - r )) • dS. (88) 

J A 

Since 



VG(r-r ) = G(r-r )(-^ + ^ ~ G(r — r )k r " 



r-r | 2 |r — r | / |r - r 



when k\r — r | 1 (that is, |r — r | A), we get finally, 

V(r ) = -i f V(r)C(r - r ) (kj—^r - k in ) • dS. (89) 

J A V l r - r 0| / 

Noting that the normal of the surface A is pointing away from the point 
r we infer that 

k- - k in ) ■ dS = k(cos6 + cos#; n ) dS, 

l r - r o| / 
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where 9 is the angle between r — r and the normal of the surface dS at 
the hole. The new thing obtained here beyond the Fresnel-Huygens principle 
used above is the geometrical factor cos 9 + cos 9- m . The Kirchoff method gives 
often good approximations though it is not entirely mathematically consis- 



tent 



19l . sect. 9.8]. Indeed, the assumption of the vanishing of ip and dip/dn 
at the screen is no longer necessarily valid for the Kirchoff solution itself. 
Sommerfeld and Rayleigh have introduced slight modifications (employing 
the mirror principle) that remove these inconsistencies. There are also a few 
exact solutions of the diffraction problem - the half plane screen and the 
wedge shaped screen - presented by Arnold Sommerfeld and others, which 
are of great theoretical interest though approximate methods will typically 
suffice in practical problems. 

5 Radiation from antennas 
5.1 The Hertz dipole 

If we consider the situation of an antenna in free space, then J ^ only 
in the region of the antenna (and the RF-circuit of course); elsewhere the 
magnetic fields satisfy the free space equations. It is useful to consider a 
simplified case where we have short piece L of a wire as an antenna at r , 



directed along the ^-direction 18]. We suppose that the current density in 



this antenna may be described as 

J(r,t) = 5(r-r Q )Ll cosut, (90) 

using the Dirac (generalized) delta- function S(. ..) in three dimensions, 
S(r) = 5(x)8(y)6(z). In denotes the constant current amplitude directed 
along the ^-direction. Inserting (l90l) into (1271) gives 

A(r 1 t) = ^° l LI ° , cosuj(t- |r-r |/c), (91) 
47T |r — To | 

from which we can calculate the magnetic field H = V x A//x/i , 

1 v / co,„( f -|,-r,|/cA x = 

4vr V |r-r | / 

cosu(t — |r — r \/c) uj sxaujit — |r — r |/c)\ , . (Lin) 

1 r ) x 



|r — r | 3 c |r — r | 2 / ' 47r 

If we use R = r — r , we can see that the cos-term in (l92l on the right 
hand side falls off as R~ 2 while the sin-term falls off only as i? -1 with growing 
R. Thus, for R~ 2 < {uj/c)R~ 1 ; that is, for (the "radiation field region") 
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R>1 (93) 



the second term (radiation field) will dominate giving 



i , uj smu(t — R/c)^ Lin . . 

H(r,t)«- L " R»J (94) 

This can be interpreted as a planar field propagating in the direction 
R = R/i? with an amplitude 

u LI sin 6 

H ° = - C ^R-> (95) 

where 6 is the angle between R and In (= the ^-direction). As the antenna 
is along the 9 = direction we can see that there is no radiation in this 
direction. 

The interpretation of our results is that the oscillating current in a short 
wire antenna generates a field which, far away form the source, approximates 
a plane wave with an amplitude which falls off with distance as 1/R. Since 
the electric and magnetic fields are related by ( ITTj) (second equation), it 
follows that the electric field amplitude E will be orthogonal to both H 
and R; in fact, we have E = i]H x R0 

Inserting (l95l) into the expression for the Poynting vector S pTT ). and 
integrating over the surface of a sphere of radius R, centered at the antenna 
(r ), we obtain, 



/ 

J\r— r h 



S-dS = -^-(^\ , (96) 



'|r-r |=B 12?r 



7 The near-field region is of interest when one considers cases of energy transfer not by 
radiation but through inductive couplings (25| . 

8 One can work out the electrical field directly (see e.g. [20, sec. II]) from the relation 

= ^ Qt 

by calculating cj) from the charge distribution g which in turn is obtained from the 
continuity equation 

Another procedure is to first determine the magnetic field from H = V x A/fi/in and 
then the electric field from E = iV x H/uieen, which follows from the Maxwell equation 
VxH = ee d~E/dt in the current-free region (J = 0). 
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for the average (over time) power radiated by the antenna. If we compare 
this to the average power RIq/2 dissipated by a resistor R through which 
flows an oscillating current Iocos(ui), then (l96l) suggests that the antenna 
can be associated with a radiation resistance 

2r] [ujL\ 2 2tt (L\ 2 /L x 2 



n °=wA-) ="tUJ am \v a - (97) 

One implication of this analysis is that, if an AC-circuit contains wires 
whose lengths approaches the order of A = c//, where / is the frequency 
of the AC-current, then one has to take into account that the wires may 
radiate a significant power and have to be treated as antennas. However, 
if the typical dimensions D are much less than A, then the system can be 
treated as an "ordinary" point-to-point circuit where Kirchoff 's laws can still 



be applied 55]. 



The above model based on (1901) . and the assumption L <C A, is referred 
to as the Hertz dipole. It already shows some general features. Thus, in the 
"far away region" the magnetic and electric field amplitudes fall off as Rr 1 , 
and since the radiated power per unit area is proportional to the square of 
the field amplitude, it falls off as R~ 2 . This is logical since the power flowing 
through the spherical surface area 4-iiR 2 must be constant (in empty space) 
and independent of R since energy is conserved. We also found that the 
direction of the electrical field (direction of polarization) lies in the plane of 
the antenna and the radius vector R; more exactly, E = E 9. 

If the antenna is of the size of the wave length, L ~ A, then the variation of 
the current along the antenna becomes an important problem. Basically one 
would have to solve the field equations for the propagation of the field along 
the antenna. The conductors act as waveguides directing the fields along 
the surface of the conductors. For good conductors there is no field inside 
the conductor and thus no energy is transported in or out the conductor. If 
there is slight resistance in the conductor then some of the energy transported 
by fields will flow into the conductor and will be dissipated as heat. In this 
sense the antenna does not as such radiate, instead it guides the fields along its 
surface. This is vivildly illustrated by the diagrams computed by Landstorfer 



and coworkers [271 ] which show how the electromagnetic energy flows along 
the surface of antennas and is spread into the surrounding space. 

If an oscillating voltage source is connected to two parallel wires it will 
generate a traveling wave between the wires which will detach from the free 
ends of the wires and radiate into the surrounding space. By bending the 
wires at the free end forming a T-dipole the part of energy radiated into 
space may be increased. The problem of determining the radiating field thus 
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becomes a boundary value problem, with excitation voltage given at the 
feedpoint, and with the boundary condition that the tangential electric field 
vanishes at the surface of the conductors. However, linear wire antennas are 
often treated by the simpler method of making some more or less well justified 
assumptions about the current distribution in the antenna, from which the 
fields are calculated as exemplified by the Hertz dipole case. The problem 
of determining the antenna current given the exciting voltage is called the 
antenna excitation problem. Mathematically it leads to integral equations 
that are numerically solved using various discretization procedures, such as 
the moment method and the method of Galerkin (l3| . For a freely available 
software (4Nec2, "numerical electromagnetic code") for computing antenna 
fields and parameters see 49, 0] - 



5.2 Dipole antennas 
5.2.1 Antenna excitation 

The dipole antenna consists of two wires extending in opposite directions 
with the feeding point at the meeting ends. A simple dipole can be made 
from a coaxial cable by exposing a A/4 end of the inner wire and make a 
sleeve of the outer conductor as shown in the adjoining figure. The added 
advantage of the sleeve is that it acts as a so called balun. In order to 
study the current distribution in a common dipole antenna we assume that 
the wires are cylinders of radius a < A, oriented along the z-direction, and 
separated by a small gap where an oscillating electrical field E z drives the 
antenna current. For a thin wire dipole- antenna it is typically assumed (based 
on the thin-wire approximation; see below) that the current distribution is 
sinusoidal along the center of the antenna (x = 0, y = 0), 

3(z) = I m sin - k\z\^j 5(x)S{y)z. (98) 

(Note that the feedpoint current amplitude I(z = 0) = Iq is related 
to the maximum current amplitude I m by I — I m ■ sm(kL/2).) Here the 
antenna extends from z = -L/2 to z = L/2, and f|98l) is consistent with 
the requirement that the current must vanish at the endpoints. (As usual 
k denotes the wave-number 2ir/X.) The approximation (|98l) for the current 
however is inconsistent with the fact that the current is concentrated on 
the surface of the conductor. Still, for calculating fields far away from the 
antenna (in terms of the radius a of the wire) the difference will be small. 
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The mathematical form of the antenna ex- 
citation problem for wire antennas is usually 
developed as follows. Differentiating the sec- 
ond equation in (1231) with respect to time, and 
using (1251 ) to eliminate (f) we obtain 



V(V • A) 



1 d 2 A 



c 



1 <9E 



(99) 



Because of the cylindrical symmetry we 
may assume that A = (0,0, A z ). Secondly we 
express A function of the current J using 
(1271) , leading to the equation 



I 



&A,(z,t) l &A,(z,t) 

dz 2 c 2 dt 2 ( ' 

?!* f L/2 f 92 1 92 ] J ( w ' *) a 

47T J_ L/2 \ dz 2 C 2 dt 2 J ^{y. - Z )2 + a 2 U 

1 dE z {z,t) 
c 2 dt 

on the cylinder (wire) surface (to repeat this is strictly speaking incon- 
sistent, see e.g. [20]). Finally, using the time-harmonic form 



I(u,t) = I{u)e- luj{t - R '^ (R = y/(u- zf + a 2 ), (101) 
E z (z,t) = E z (z)e-^, (102) 

we end up with Pocklington equation (which does not appear in this form 
in H) 



I(u) l-^ + k 2 } =du = (103) 

! \dz 2 J An^/(u - z) 2 + a 2 V ' 



L/2 



-ieeouEJz) 



where we have used ( fl"4l ) and the fact that u/c = k. This is an integral 
equation in the unknown current / in terms of the excitation field E z . Since 
it solves for the current /, given the input potential V, it gives the impedance 
Z = V 1 1 of the antenna. If the antenna is assumed to be a perfect conductor 
then E z = at the surface of the antenna since a non-zero tangential field 
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would generate an infinite current (more precisely, in a perfect conductor it 
generates a field that cancels it such that the total tangential field becomes 
zero - more on this in the section on shielding). E z differs from zero only at 
the feedgap where it may be assumed to be —V/A, where V is the voltage 
fed into the antenna, and A is the gap size. The same equation applies 
also to the receiving antenna in case which E z = —E^ (in the right hand 
side of Eq. (ll03l) ) at the surface of the conductor; i.e., the scattered field 
offsets the incoming field such that the total tangential component vanishes. 
Pocklington developed a somewhat more general theory of wire antennas 
in |38j as he considered wire antennas of arbitrary shapes, such as circular 



rings and helixes. His basic idea was to consider the wire antenna as being 
composed of an array of infinitesimal Hertz dipoles. The total electric field 
can then be obtained by summing (or integrating) the contributions from the 
elementary Hertz dipoles, each characterized by a current J(s) and direction 
(is, 



E(r) = -V J dGir I(u)du + k 2 J I(u)G(y - s(u))^du, (104) 
where 

G(r) = wo^-. (105) 

The integration in Eq. (1104l) . which can be derived from Eq. (11440 and 
Eq. (11451) . is along the wire antenna {u is the length parameter along the 
curve defined by du = \ds\). Finally Pocklington imposed the boundary 
condition that the tangential component of the electric field ( 11041) vanishes 
at the surface of the conductor. In the limit of a vanishing small radius this 
leads to the equation^! 



9 For closed antenna loops, or where current is assumed to be zero at the ends of the 
antenna, the first term in (|104|) becomes by partial integration 



V / G(r- S ( M ))^<i«. 
J du 



Pocklington then evaluates l|104p at surface of the conductor of a small radius 5. For 
small S the leading contribution to the field at a point on the surface is from the nearby 
section of the conductor yielding a term proportional to ln(M/<5) (for some constant M). 
Indeed, note that dzj\J z 2 + S 2 — > ln(2C/<5) for small S. The electric field parallel to 
the conductor at its surface thus becomes proportional to 



d 2 I(u) 
du 2 



+ k 2 I{u) \n{M/5) 
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^ + k*I(u) = 0, (106) 

i.e., the current is a sinusoidal function of the wire length parameter 
u. (This equation follows from the approximation of ([271) . that A z on the 
conductor is proportional to the current at that point, in combination with 
Eq. (1107p . For a pedago gic discussion of various approximation schemes for 
the antenna current see [30|.) 

If we return to the version fll03f) and let A — > we arrive at the deltagap 
case where the feedgap is assumed to be infinitesimally small at z = 0. Spe- 
cializing the earlier considerations of the dipole antenna to this case (usually 
associated with the name Hallen [15]) we get the equation 

^ + k 2 A z = (z^O), (107) 
which suggests that A z is of the form 

A z = -(Bcos(kz) +Csm(k\z\)), (108) 

where B and C are constants to be determined. This form ensures the 
symmetry A z (z) = A z {— z) and I(z) = I{—z). When differentiating the 
| z | -argument this introduces a discontinuity at z — 0. In fact, using the 
equation 

- 1 E - = i£ + k ' A " {m 

in order to compute — V = J E z dz we obtain the value 2Ci which deter- 
mines C in terms of V. B finally is determined by the condition that the 
current J, solved (numerically) from 



A Z = -(B cos 

should vanish at the ends of the antenna, I(±L/2) = 0. Having solved 
for I(z) we obtain the impedance Z of the antenna from Z = V/I(0). (If we 
use the approximation that A z is proportional to the current I at the wire 

when neglecting terms small compared with ln(M/<5) as S — > 0. The vanishing of the 
tangential field component therefore leads in the limit of a thin wire to (| 106|) excepting the 
end points and feedpoints of the antenna. It is to be emphasized that the above argument 
does not demonstrate that sinusoidal currents do produce zero tangential fields for nonzero 
radius. In fact, for sinusoidal current the tangent field will in general look sinusoidal too. 



(kz) +Cwn(k\z\)) = /ifi / I{z')-^dz' 



(110) 
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surface, then the condition 7(±L/2) = would simply become B cos(kL/2) + 
C sin(/cL/2) = 0, but such an approximation is too crude in order to yield a 
reliable estimate of the impedance of the antenna.) 

In the integral in Eq. (11101 ) it is assumed that the current I is at the axis of 
the antenna leading to the simplified "kernel" e lkR /4nR used above (called the 
reduced kernel approximation) , with R denoting the distance from a surface 
point to a point on the axis (see Eq. (llOip ). In the cylindrical antenna theory 
the reduced kernel 



This kernel is based on the replacing an axial current / at the cross section 
z with a surface current I/2ira uniformly distributed on the circumference at 
the section z. One sees that (11101) is, in fact, mathematically inconsistent if 
one uses the reduced kernel, in the sense that the left hand side is continuously 
differentiable at z = (supposing the current I is "well behaved") while the 
right hand side is not [13fl. The cylinder kernel fares better on this point since 
it has a singularity R = (for the reduced kernel one has always R > a). 
Indeed, the kernel K(z) grows like — ln(z) when z goes to zero. However, 
integrals J K{z — u)I{u)du still remain finite because f\n(z)dz = z In z — z 
stays finite when z approaches 0. 

Returning to Eq. (1 11 Oft we can discretize it by dividing the antenna into 
N = 2Q + 1 slices of length A = L/N with the feedgap at z = (corre- 
sponding to the index n = Q). The midpoints of the slices are given by 
z n = (n — Q)A. A discretized version of (11101) can thus be written as the 
matrix equation 



gifc-y/ (u— z) 2 +a 2 



47Ta/(u - z) 2 + a? 



in Eq. (I103R becomes replaced with 




(111) 



N-l 




(112) 



k=0 



where 
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pn pA/2 ik^/((n-k)A-z) 2 +4a 2 sin 2 § 

M n , k = / / . dycfe, (us) 

4vr Jo y-A/2 J(( n _ fc)A - z) 2 + 4a 2 sin 2 f 

Vfc = ^(Scos(A;^) + Csiii(A;|^|)), (114) 
C = ilA/2 (U is the applied feedgap voltage). (H5) 

The matrix equation V = MI is easily implemented with mathematics 
software such as Mathcad, Matlab, etc. One solves for I = M~ l V and 
adjusts the coefficient B such that 1$ = Ijsr-i = 0. Writing the vector V 
in (dUl) as V = (B/c)V^ + {C/c)V® one obtains B from B{M~ X V^) Q + 
C(M~ 1 V( 2 >) = (which is thus generally a complex quantity). The ad- 
joining figure (FigH|) shows the results for imaginary and real parts of the 
current in the case Q = 100, a = 0.001A and L = 0.48A (A = 12.5 cm). 
As can be seen the calculated current deviates slightly from the sinusoidal 
current. The impedance was obtained as Z = U/I{z = 0) = 73.4 — i4.3Q 
whose real part may be compared to the radiation resistance (73 fl) for the 
half- wave dipole to be calculated below assuming a sinusoidal current. This 
result indicates that by making the dipole a bit shorter than A/2 one can 
make the complex part (reactance) of the input impedance disappear, which 
simplifies the impedance matching of the antenna (such an antenna is called 
a "resonant antenna"). 
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Fig. 1: The figure shows the distribution of the current along a dipole an- 
tenna, with an applied deltagap potential of 1 V at the center. The 
current has been calculated by solving the discretized Hallen equation 
(jll2l) . The solid line represents the real part of the current, the dashed 
line the imaginary part, while the dotted line represents a pure sinu- 
soidal current. Parameters: L = 0.48A, antenna radius a = 0.001A; 
the calculated input impedance becomes Z = 73.4 — i4.3 O. Note that 
in the literature it is often the absolute magnitude of the current that 
is display instead of the real and imaginary part separately. 



Another approach to solving the Hallen equation (11121) is to express the 
current I as trigonometric sum 32 



I(z) = J2D n sm — (-- 

n=l ^ ^ 



(116) 



The idea is thus to try to approximate the current using basis functions 
that are expected to be similar to the true current. We note that (11161) 
automatically satisfies the condition 7(±L/2) = 0. If the expression ( 11161 ) is 
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inserted into Eg. (II 101) we obtain an equation of the form 



J D n sin ^— ^— — \u\\ j K(z — u)du = B cos{kz) + C sm{k\z\). 

(117) 

Choosing N + 1 points in the interval (-L/2, L/2) we can write (I117p as 
a matrix equation 



N 



J2 K m,nD n + K mtN+1 B = Csm(k\z m \), (118) 



n=l 



where 



K m ^ n = ^ J ^K(z m - u) sin (^(^-\u\^du (n < N), (119) 

K m:N+1 = -COs(kz m ). (120) 

As before C is related to the applied feedgap voltage U by C = iU/2, 
and K(z) denotes the kernel (jllip . When selecting the points z n one should 
take care not to make the (N + 1) x (TV + 1) matrix K singular. (This may 
happen if the points are selected symmetrically along the 2-axis with respect 
to z = 0.) Using N + 1 = 6 points we obtain a five term trigonometrical 
approximation of the current which hardly improves much going to higher 
values of N - the computational cost increases rapidly due to the integration 
over the interval (— L/2, L/2). 

5.2.2 Receiving dipole antenna 

For simplicity assume that the incoming field is parallel with the dipole 
antenna, chosen again to be along the z-axis. For a perfect conductor the 
incoming field will induce a current and a reflecting field E such that the 
tangential component of the total field E + E^ is zero. From this we obtain 
that E z = —E z m \ For an incoming plane wave we can assume that is 
constant along the antenna0 We can therefore write the solution to Eq. (11091) 

cLS, 



10 In the more general case the last term in (|12ip is replaced by 

- r^")(C)sin(fc(z-C)R- 
^ Jo 
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A z = -(Bcos(kz) + Csm(k\z\)) + -E ( / n) (A/2 < \z\ < L/2). (121) 

C UJ 

Making use of the Lorenz gauge condition 

V A + ^ ^ ^ 2 ^0 

c 2 dt dz c 2 

the potential can be expressed in terms of A z , 

, c 2 dA z 
^ = " — ■ 

For a thin wire dipole antenna we can use the "Pocklington approxima- 
tion" and assume that A z is proportional to the current [3J|. In the open 
circuit case we have the condition that the current is zero at \z\ = A/2, L/2. 
The consequent condition A z {z) = for \z\ = A/2, L/2, then leads in the 
limit A = to the solution, 

B = —Ef n \ c= ^ ^-^kL,2) ^ 

uj uj sm(kL/2) v ' 

In the same limit A — > we then get for the induced open circuit voltage, 



if a / \ c 2 fdA z (z = A/2) dA z (z = -A/2)\ , . 

u=*w*)-<K-±m = to\ \z — dz ) ^ 

_2E { z in) l-cos(JfeL/2) 
~~ ~~k sin(A;L/2) 

Especially we obtain for the half- wave dipole (L = A/2): U = EzX/ir. 
5.2.3 Dipole field 

In order to be able to derive some analytical results for the dipole antenna 
of length L = 21 we will assume a sinusoidal current (l98l) (which, as pointed 
out above, is typically sufficiently accurate for estimating the far field.). As 
in the case of the Hertz dipole we can then calculate the radiation field. 
For an antenna oriented in the z-direction the electric component of the 
radiation field (in complex representation) in the far-field region (r > A, r is 
the distance from the antenna and A is the wavelength) thus becomes 

= § cos(klcose)-co S (kl) 

sin Aur 
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Here Iq is the amplitude of the sinusoidal antenna current, 9 is the unit 
vector in the direction of the ^-rotation, k = 2n/X is the wave-number, and 
r] is the "wave impedance". In Eg. (11241 ) we have ignored the effects of the 
surrounding, such as the circuit board (PCB) of the transmitter. In our 
measurements it was found that the PCB contributed to an anisotropic form 
of the radiation field in the horizontal plane which can be cured enclosing 
the PCB in a symmetrical (cylindrical) shielding box. Eq. (11241) is a special 
case of the more general form 

E(r) = ^-G(M) (125) 

of the far-field electric component. Since the power flux of the radiation 
field is proportional to E 2 , Eqs. (11241 1125H predict that the power depends on 
the distance r as r~ 2 . Power P is usually measured in the units of decibel 
(dB) relative to a standard power P (such as 100 mW), 

P dB = lO-log lo 0^. (126) 
which is equivalent to the field strength. 

Thus, if we have a r~ 2 -dependence of the power then it is reduced by 
20 dB for every decade of distance and by ca 6 dB for every doubling of 
distance. From this relation one could deduce that, if the power at r = 4 m 
is -10 dB, and we measure a power -40 dB at an unknown distance r, then 
r must be 4 m x lo^ 40-10 )/ 20 £s 126.5 m. Besides the distance r the received 
signal strength depends also on the antenna orientations. 

A common dipole antenna is the half-wave 
(A/2) T-dipole with I = A/4 and kl = tt/2. 
In this case the field ( 11251 ) component Eg be- 
comes, 

E e = ^ COS ^ COS V. (127) - 

2-rt sin 9 — 

The total radiated power is obtained by in- 
tegrating Ej/2ri over a spherical surface r = 
constant, 

1 fvImV f w/2 f COS U COS 9)Y 



/ 



5 Radiation from antennas 



49 



The integral cannot be solved analytically but can be evaluated numeri- 
cally on the computer. The end result is that, 

P tot « 0.097 • rj ■ I 2 m} (128) 

and identifying this with R s Iq/2 (note that for an A/2-dipole I m = Jo) we 
obtain for the radiation resistance R s of the A/2-dipole, R s ~ 2 ■ 0.097 • rj ~ 
73 Q. We can generalize these calculations using (11271) to antennas of lengths 
L n = (2n — l)A/2 (the calculations are simplified for these special lengths) 
which yield for the corresponding antenna resistances the expression 



P n = 30fi-Cin((4n-2; 



7T 



(129) 



where the cosine integral is defined by 



Cin(x) 



cosy 



y 



dy. 



(130) 



A variant of the dipole antenna is the folded A/2-dipole where the end- 
points of the T-dipole are connected by a wire. Thus the current will counted 
twice when evaluating the radiation field which therefore will be twice com- 
pared to the T-dipole. The radiated power will consequently be four-fold. 
Hence the radiation resistance of the folded dipole is also four-fold, or about 
4 x 73 Vt = 292 a 

In a bit more careful treatment of the folded 
dipole case 0, section 9.5] the conventional 
analysis is based on the trick of viewing the 
system as consisting of two parallel coupled 
dipoles separated by a small distance d (see 
figure). These are excited by voltages U/2 
and U/2 (symmetrical case (a)), and U/2 and 
—U/2 (asymmetrical case (b)). The current 
I of the original problem is then obtained as 
the sum of the currents I a and h (at the feed- 
points) for the separate cases. In case (a) we 
have the equation 



u 

2 



U 



U 
2 



(b) 



u 

2 



U/2 = Z u I a + Z 12 I a , 



and assuming that for small separation d the mutual impedance Z 12 is 
close to the self-impedance Z = Z\\ of the single dipole we obtain, 
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I, 



a 



u 

~AZ 



The asymmetrical case (b) can be handled with the help of the transmis- 
sion line theory (see Appendix IDI). The input impedance of two- wire cable 
of length L/2 short circuited at the end is Z b = —iZ c tan(kL/2) where Z c is 
cable impedance. Thus the current lb becomes 



When L = A/2 we obtain that the folded dipole impedance Zf is four 
times the simple dipole impedance, Zf = 4Z. In this case only I a contributes 
to the total current since the short circuited A/4-length transmission line has 
infinite impedance which suppresses the asymmetrical mode (b). 

In the cases considered we have assumed that the feedpoint of the dipole 
antennas is centered at the midpoint (corresponding to z = 0). If the 
feedpoint is off-center by an amount h, then the current at the feedpoint 
- assuming a sinusoidal current I m sm(kL/2 — k\z\) for the antenna - will 
be I m sm(kL/2 — k\h\) which for a A/2-dipole becomes Ih = I m cos(kh). 
The radiated power can therefore be written P = 1/2 • 3ft (Zi^) = 1/2 • 
9ft(Z/ cos{kh) 2 )l1 implying that the off-center input impedance is Zh = Zj cos(kh) 2 . 

5.3 Mutual impedance 

A current in a conductor 1 generates an electric field which induces a current 
and its reaction field in any nearby conductor. This mutual influence is the 
basis of the use of antennas. The mutual influence is measured in terms of 
the mutual impedance. If we consider two (unloaded) dipole antennas and 
denote the voltages and the currents at the feedgaps (antenna terminals) as 
Vi,Ii, then we have the relations, 



lb = — -. — iZ c t&n(kL/2). 
The impedance of the Zf of the folded dipole is therefore given by 
1 I a + I b 1 i 



Z f U AZ 2Z c tan(fcL/2)' 



Vi — + Z12I2, 
V2 = Z % \I\ + Z22I21 



(131) 
(132) 



where Z n ,Z 22 are the antenna impedances, and Z 12 = Z 2 \ is the mu- 
tual impedance (see Appendix and the Reciprocity Theorem). Following 
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Fig. 2: Field by antenna 2 induces an open-circuit voltage V\ over the termi- 
nals of antenna 1. 

Bechmann [J] one can compute the impedance for a set of N conductors by 
starting from the Poynting's energy theorem (fl~9lh Taking the time aver- 
age for harmonic fields the term d/dt{- ■ ■ ) disappears, and we obtain using 
phasors, 

Pav = -\J^(V-y)dV. (133) 

Denote by E» the electric field generated by the current density Jj of the 
zth conductor. Then the equation (11331) becomes 

p ™ = -\Y.J^ ( Ei • J i) dv - ( 134 ) 

The term / (E< ■ 3* + Ej ■ J*) dV measures the mutual influence between 
conductors i and j. Since Ej can be supposed to scale with the amplitude 
of the current density Jj, due to the linearity of the Maxwell equations, one 
can define a current-independent quantity called (mutual) impedance by 

Zl3 = 2ljJ( El ~ + Ej ' J ^ dV ' (135) 

where Ii denotes the characteristic current amplitude of the ith conductor 
(one may refer to the current amplitude at the antenna terminals, or the 
maximum current amplitude of the conductor). 
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In the Appendix on the Reciprocity relations we consider the transmission 
between two antennas. If the antenna 2 is a transmitting dipole while antenna 
1 is the receiving dipole then it was demonstrated that the open circuit 
voltage VI induced by the field generated by 2 is given by V\ = Z12I2 (see 
Eq.Q215J) in the Appendix) 

V 1 = jJe 2 - JidV. (136) 

This is typically calculated by inserting a sinusoidal current for Ji, letting 
E 2 be the incident field and integrating over the antenna. 

In the situation considered by Bechmann E2 would refer to the incoming 
field and Ji to the induced current, and the total field (satisfying the BC) 
would be E x + E 2 . Then the integral in Eq. (11361) over the antenna can be 
interpreted as an electromotive force generated in the antenna0 This voltage 
can furthermore be expressed as V\ = hi (9) ■ E 2 where hi(9) defined by (for 
a linear antenna) 

1 r L/2 

hi(9) = — / e ikzcos9 Ii(z)dz (137) 

h J -L/2 

is called the effective antenna length (9 is the angle between the wave 
direction and the linear antenna - see illustration). The exponential fac- 
tor comes from the fact that the far field E 2 (r) varies as exp(ik • r) = 
exp(ik • rn) exp(ikz cos 9) along the receiving antenna. Again using sinu- 
soidal currents as an approximation for a A/2-dipole we obtain 

11 Carter Q argues that the impressed field E z {z) at z of the antenna induces a voltage 
E z (z)dz at the point z over the section dz, and a current I(z), and that E z (z)dz I(z) 
by reciprocity is equal to I(0)dV where dV is the corresponding induced voltage at the 
antenna terminal z = 0. Summing over all sections gives the induced total voltage V(0). 

In [48] the authors adopt the approach to let E 2 be an arbitrarily prescribed incoming 
field which generates a current (7E2 in the antenna which is assumed to have finite con- 
ductivity a. This current is then added to the Maxwell equations for the "reaction" fields 
Ei, Hi in the antenna (conductor): 

V x Ei + i/i^owHi = 0, 

V x Hi - crEi = crE 2 . 

(In the second equation the displacement term has been dropped as negligible for a 
good conductor.) Combining this with the equations for the fields outside the conductor 
(free space or a dielectric) leads to a system of equations describing waves traveling along 
conductors which was treated by Sommerfeld already in 1899 (a more recent reference is 

S3). 
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[■l/2 2 r /2 

h{9) = e lkzcose cos(kz)dz = - cos (u cos 9) cos udu (138) 

J -L/2 k J 

2 cos (| cos 9) 
k (sin#) 2 

Hence for a A/2-dipole at 9 = n/2 we obtain h = \/ir. This agrees 
with the result in section 15.2.21 As discussed in an earlier example, a near 
planar EM-field with power flux of 100 mW/m 2 corresponds to an electric 
field amplitude of E = 8.7 V/m. If this is aligned with a receiving A/2- 
antenna (i.e. the antenna is parallel with the polarization), and given A = 
0.125 m, the field will induce an open circuit voltage of amplitude h x 8.7 
V/m = 0.125/7T x 8.7 V = 0.35 V between the antenna terminals. 



5.4 Antenna above ground 

We have discussed the antenna radiation field in the open space. However, 
antennas are usually situated near the ground which will thus affect the radi- 
ation field. As shown in section [3761 the static field of a charge in presence of a 
dielectric or conducting plane can be handled using the image method; that 
is, the total field above the ground is the sum of the direct field and the one 
generated by an opposite "mirror" charge. An analogous situation prevails 
also in the case of moving charges in antennas. Thus, if the transmitter- 
antenna T is situated at r = 0, and the ground is on the level z = —h, then 
the radiation field is the same as in the free space case with a second trans- 
mitter T M situated at r = (0,0, —2h). T M is a mirror image of T reflected 
in the plane z = —h and is fed with the current Im = —I- For a verti- 
cally oriented dipole antenna above the ground the resulting far field E tot (r) 
becomes, 

E tot (r) = E(r) + E(r-D), (139) 

where D = (0, 0, —2h) and E is the dipole field of Eq.(J1240. This equation 
for the total field is justified in the case of a perfectly conducting ground since 
it satisfy the condition of a zero tangential total field at the interface. 
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The field E tot represents the sum of the direct field and the reflected field. 
One thus speaks of a "two ray model". If we use Eg. (11391) and Eg. (11241) in 
order to calculate the power (which is set proportional to Ef) as function of 
the distance from the transmitter we obtain Figj3l The solid curve represents 
the case when the receiver and the transmitter are at the height h = 1 m 
above the ground. As can be seen the interference has a marked effect and 
the power does not decrease monotonically with distance as it would do in 
the free space. In contrast, the dashed curve representing the case h = 0.1 
m shows a monotonic fading with distance. 

Our measurements in an outdoor soccer field revealed that the ground 
(sand, gravel) was best described as a dielectric with the relative permittivity 
e ~ 3. In this case we need to employ the Fresnel relations for the reflected 
part of the horizontally (H) and vertically (V) polarized components of the 



field, ( 1581) . (1591 ). Eq. (11391) has then to be modified by adding the reflexion 
factor, 

E tot (r) = E(r) + p„(0)E v (r - D) + p h (6)E h (r - D), (140) 

for the horizontally and vertically polarized components. (For vertical 
antennas E v = E .) Eq. (11401) can be justified by the boundary conditions at 
the interface when the field can be assumed to be approximately planar. A. 
Sommerfeld has obtained an exact (fo r the idealized model) solution for the 



Hertz dipole above the ground [43|, [4j|. An alternative (but mathematically 



equivalent) approach has been presented by Weyl [53|] where he decomposes 
the dipole field into a sum (integral) of planar waves. Thus one may apply the 
Fresnel relations to each planar wave component an then sum the components 
which yield the solution. 
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Fig. 3: Variation of power with distance due to interference calculated accord- 
ing to the model of Eq. (ll39l) . Antennas (transmitter and receiver) are 
supposed to be vertically oriented. Solid line corresponds to antennas 
at the height h = 1 m; dashed line to h = 0.1 m; dotted line to /i sen der 



0.5 m and h r 



0.1 m. 



5.4.1 Sommerfeld's analysis of the antenna above ground 

We consider a vertical Hertz dipole at a height h above a planar ground 
characterized by parameters e and a. Above ground the radiation field can 
be considered as a sum of the primary (open space solution) excitation plus a 
reflected part, while the field in the ground consists of the transmitted part. 
The task is to is to determine the reflected and transmitted fields by invoking 
the boundary conditions at the interface. The difference here compared to 
the earlier reflection problems is that we no longer assume that the primary 
field is a planar wave but is given by the full Hertz dipole field (in complex 
notion), 
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e iuj(r/c-t) 

A(r, t) = fifiollo— , (141) 

where r = |r — ro| is the distance from the antenna. The electric field is 
then given by E = — V0 — dA/dt. However, the 0-part can be avoided if we 
use the Maxwell equations and the Lorenz gauge to derive the equation 

A s 1 d 2 A IdE , . 

V(V-A)-^ = ? ^. (142) 

Then for harmonious waves we can replace d/dt(- ■ ■ ) by — ■ •), and 
introducing the Hertz function n through 

A = 317, (143) 



we obtain 



c 2 dt 



V(V • n) + k 2 U = E. (144) 



(As usual kc — u.) In open space we have the Hertz dipole solution which 
we write simply as (the important part here is only how it depends on r) 



n = z , (145) 

r 

for a dipole oriented along the vertical direction z. It satisfies the free 
space wave equation 

V 2 II + k 2 fi = 0. (146) 

It will be assumed that the dipole is .at z — h while the surface of the 
ground corresponds to z = (and air for z > 0). Eq. (11451) represents the 
primary field n prim excited by the antenna. With no ground this would 
constitute the total field. However, with the ground present the total field 
will become a sum 

n = n prim + n S ec, (147) 

where n sec represent the secondary "reaction" field generated by interac- 
tion with the ground. In case of the perfectly conducting ground we already 
found the solution for n sec through the image charge principle. Now we as- 
sume that the ground has finite conductance a and an index of refraction 
given by n% — e + ia/ueo. Earlier in section [3751 we discussed reflexion and 
transmission for planar fields and how the reflected/transmitted part was 
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found by considering the boundary conditions of the field at the interfaces. 
The same applies here. For this end one may take advantage of the cylin- 
drical symmetry of the situation and employ Bessel functions. In terms of 
cylindrical coordinates (p, <p, z) the distance r from the antenna is given by 

r 2 = p 2 + {z _ ^2 

From the cylindrical symmetry it follows that n prim and Il sec have com- 
ponents only along the z-direction for which we also write n prim and II sec . 
Now the primary field can be expressed in terms of the Bessel function J as 



an integral |44j, p. 243] (known as Sommerf eld's equation) 

Akr poo 



ikr poo j- if 

n prim (p, z) = — = Jo(pC) e-^-M (148) 
r Jo 7 

for the region z > 0, and where 



(The real part of 7 is taken to be positive in order to ensure the conver- 
gence of the integral.) This suggests the ansatz 



roc 

U 8ec (p,z) = F(()J (pOe-^ +h U(, (149) 
Jo 

POO 

II earth (p,z) = / F E (()J (pC)e^ h d(, (150) 
Jo 

for the secondary field and the transmitted field in the ground. Here 



1e = VC 2 - k%, 

with k E = n E k (the term jh in the exponents is added only for conve- 
nience in dealing with the boundary conditions). n ea rth represents the field 
transmitted in the ground. The functions F and F E will be determined by 
the boundary conditions for the electric and magnetic field at the interface 
z = 0. The conditions are that the tangential components E p and H v shall 
be continuous across the boundary z — 0. These components are given in 
terms of II by, 



E > = TzW (151) 

H, = (152) 
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and similarly for the earth fields (z < 0). The boundary conditions at 
z = can thus be written, 



d_dU_ _ d dn carth 
dz dp dz dp 

]? 2 ^ _ ,2 dearth 

(We have assumed that the magnetic permeability is p = 1 also for earth.) 
The first condition can be integrated over p leading to 

<9n _ dn carth 

dz dz 

Inserting the expressions for II = n prim + Il sec and n ear th one is finally led 
to the following equations for the functions F and F E , 



F — — ( 1 — 2 2lE \ , (155) 

F E = 2 2( \ . (156) 

n E l + 1e 

These functions are related to the Fresnel coefficients p and r Eq. flo8l ). 
although the relation is not that transparent due to the integral form in the 
present case@ However, the solution for z > can be written (H, p. 251], 

n = — + ^ - 2 T J ((p)e-^ ^. (157) 

The second term corresponds to the reflected field in the case of a perfectly 
conducting ground where 



r* = a/p 2 + (z + h) 2 . 

Indeed, we see that the integral in Eq. (11571) vanishes in the limit of \n E \ — > 
oo which corresponds to the case of a perfectly conducting ground. Hence the 
integral term describes the deviation from the case of a perfectly conducting 
ground. 



12 As pointed out above there is an approach by H Weyl (53| from 1919 (described in 
471 . sec 9.29]) that may be more transparent in this regard, since it develops exp(ikr)/r 



in terms of plane waves starting from the identity exp(ikr) / ikr — f. exp(ikrr])dr]. 
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5.5 Breaking point, non-smooth surface 

In the case of a conductive medium the permittivity becomes a complex 
quantity, e = e' + ie" . For a typical ground the imaginary part is negligible 
at the frequencies considered here (2.4 GHz). Assuming a real permittivity 
there is a special grazing angle, the Brewster angle 9b, as explained above, 
at which the reflected vertical component disappears. If the transmitter and 
the receiver are at the heights h\ and h 2 above the ground, then the Brewster 
angle corresponds to a distance 

D B = = (h + h 2 ) V~e. (158) 

taxi 9 B 

This distance is of interest since for r > D B the reflexion coefficient for the 
vertical (V) component turns negative and approaches the value -1. (This is 
in contrast to a perfectly conducting ground where the V-reflexion coefficient 
is always 1.) From this it follows that, for vertical antennas at large distances, 
r ^> D B and r ^> 2nhihi/X, the power varies as 

(hih 2 ) 2 , , 

Poc^^, (159) 

where h\ and h 2 are the heights of the antennas above the ground. The 
"breaking point" where the r~ 4 -dependence starts to dominate is defined as 

r B p = — ^— . (160) 

Interesting it seems that the early investigators of radio communication, 
such as Sommerfeld, did not pay much attention to the drastic interference 
effect of the dielectric ground resulting in the r~ 4 power-distance relation. 

The "mirror" (or two ray) model discussed above assumes a "smooth" 
ground. If we measure roughness s of the surface as the typical height varia- 
tion of the surface, then the Rayleigh criterion for smoothness requires that 

16smt> 

for the grazing angle 9. This corresponds to the idea that such a height 
variation causes a phase difference less than 2s sin 9 < A/8; that is, less than 
45°. Along similar lines of thought we can obtain an estimate how much 
a rough surface affects the reflection in terms of a scattering loss factor f s . 
Suppose height variation s around an average height s = is Gaussian dis- 
tributed with standard deviation S s . Consider the parallel reflected rays with 
the grazing angle 9. Due to the height variation we get an extra phase factor 



5 Radiation from antennas 



60 



exp(ik2s sin 6). Summing over all the rays in the 6 direction the contribu- 
tions have to be weighed by the Gaussian distribution, which leads to the 
integral, 



f s = — — / exp I ik2ssia9 — —p^] ds = exp ^ 



(162) 

The corrected reflexion coefficient is thus obtained as / s p. Thus in the 
case the distance between the transmitter and the receiver is 10 m, the devices 
are at the height 1 m (corresponding to 6 ~ 11.3°) and 8 S = 1 cm we get f s 
= 0.98, and when S s = 5 cm we get f s = 0.62 (@ 2.4 GHz). 



5.6 Ground wave 

The two ray model seems to account rather well for the ground effects when 
comparing to measurements. It is based on calculating the reflection of rays, 
or planar waves, from the ground which usually are good approximations in 
the far field region. The mirror method gives the exact solution in case of a 
perfectly conducting ground but for a general dielectric ground the situation 
is different. Near the ground the propagation may be altered since, due to 
the air-ground interface, the EM-field might no longer be a pure transversal 
TEM field. This is in analogy to the waveguide case discussed earlier. The 
so called ground or surface wave (Oberflachenwelle) may become of interest 
in settings where the transceivers are directly on the ground. Sommerfeld in 
fact originally got interested in the problem of the field of the dipole near 
the ground in order to find out whether the dipole solution contained the 
surface wave conceived earlier by Zenneck [54]]. It may be however in the 
interest of clarity to treat the simpler case of Zenneck. We consider a plane 
wave that travels in the x-direction over a plane ground defined by z = 0. 
Thus we have a ground for z < and e.g. air for z > with the interface 
at z = 0. We assume that electrical field is polarized along z. However from 
the discussion of the skin effect we already learned that the field becomes 
damped in a conducting media and that E z therefore also is a function of z 
(besides of x), E z (z). Then from the V-E = dE x /dx + dE y /dy + dE z /dz = 
it follows that E x and E y cannot both vanish because then we would have 
dE z /dz = 0. We may thus assume that E is of the form (E x , 0, E z ) and H is 
of the form (0, H y , 0). From Maxwell equations we obtain the wave equations 
in the harmonic case, 
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d 2 E z , 2 7-1 T-i 9 E x . 

-^Y + k E z + ian/j ujE z = q-^, (163) 

d 2 E x , o t-i t-i d 2 E z , 

-Q-p + k E x + ighhquE x = q-q-, (164) 

where fc 2 = /z/ioeeou; 2 . For air we set e = // = 1 and a = 0, and for the 
ground e = e E , f-i = //# = 1 and cr = cr^. From Eq. (11641) we can again infer 
that dE z /dz ^ leads to a non vanishing ^-component. For these linear 
equations we have solutions of the form (due to their couplings E x and E z 
must have the same exponential factor) 



E x = Aexp(ik x x + ik z z), 
E z = B exp(ik x x + ik z z), 

which inserted into the equations (11631 1164D yield the relations (and their 
duplicates for ground) 

k 2 + k 2 = k 2 = k 2 + iafxfxouj, (165) 
Ak x + Bk z = 0. 

From the boundary conditions (continuity of the tangential component 
E x , and the normal component D z ) at the air-ground interface z = we 
further obtain, 



k x E ^X! 

A = A E , 
B = i E B E , 



with 



e E = e E + i • 

The index E refers to earth/ground. Summing up we have, after setting 
A — 1, altogether 7 equations for the 7 unknowns (k x , k y , k xE , k yE , A E , B, B E ) 
From these relations we obtain, 
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k 7 = ±k- 



k Z E = £e ■ k z . 



This is nothing but a solution of the Fresnel equations in case of zero 
reflexion (corresponding to the Brewster angle). Zero reflexion was here 
enforced by the initial assumption Eq. (1165l) . The k x and k z ,k Z E are in 
general complex. In order that the field stay finite as z — > ±00 we have 
to choose the negative sign of the square root in the above equation. If 
^(cb) <S 1 then exp(ik x x) results in a damping factor that can be written 

as exp ^— x ■ ky^T/2cr/ (2e uj) J which vanishes as x — > +00. These surface 

waves were studied by Cohn (1900), Uller (1903) and Zenneck (1907) 0, 
sec. 9.10]. Sommerfeld also derived a surface wave from his vertical dipole 
solution discussed above in the limit k x r — > 00. In this limit he obtains 



(z > 0) for the Hertz function a term 4% p. 256] 



n = A. e 



(166) 



where k x and fc 2 are as above (A is a constant). This describes a sort 
of a two-dimensional wave in the z-plane since it is proportional to 
instead of 1/r as in the case of "space waves". Interestingly there has been an 
ongoing controversy about the significance and reality of the surface waves 



that has lasted for over a century |5J, |43j, |53j, |33|, |34], |22|, [5jJ, QJJ. One of the 



famous "legends" promulgated as recently as 1998 in a review of the field by 
Wait (5l| is that Sommerfeld made a "sign error" in his pioneering paper [43| 
with serious consequences for the interpretation of the results. Sommerfeld 
never admitted to any such "errors" (see for instance [44j section 32 and 
ch. 23 in [3]), and the recent study H] by Collins indeed reaffirms "that 
the famous sign error is a myth". Collin traces the myth back to a short 
paper |33j] by K A Norton in 1935 which asserts that there is a sign error in 
Sommerfeld's 1909 paper - this allegation was later uncritically repeated by 
numerous authors (for instance in the book [47J) although the exact location 
of the "error" was never revealed. More importantly Norton did endorse the 
concept of the surface waves based on his "corrected" version of the theory. 
Nevertheless, the asymptotic evaluation of the Sommerfeld solution remains 
a somewhat tricky business. In fact, Kahan and Eckhart have presented a 
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series of studies [22| where they argue that no surface waves are contained 
in the dipole solutions, contrary to the opinions of Sommerfeld, Norton, and 
others. They asser10 that a careful evaluation of the Sommerfeld solutions 
shows that the surface term will be canceled by an equal term of opposite 
sign missed by Sommerfeld. Kahan and Eckhart thus agree with the 1919 
analysis by Weyl. The interesting thing is that Sommerfeld's result agrees 
with Weyl's result 0, p.937]. 

The starting point in Sommerfeld's approach (which no one diagrees with) 
is to replace Jo(Cp) m Eq. (11571) by 



MCp) = \ (h^(c p ) + tff(Cp)) = ~ (itf } (Cp) - ^ 1} (-Cp)) • 

This makes it possible to rewrite the integral in Eq. (11571) as an integral 
extending from infinity to infinity, 

n= / H£\(p)e-^ f E(dC . (167) 
Jw n 2 El + lE 

The path W extends from — oc+i5 to 0+i<5, and from 0— iS to oo—iS where 
5 > is an infinitesimally small real quantity. This rule is enforced in order to 
ensure the correct evaluations of the square roots involved in the integrand. 
The tricky part is the asymptotic evaluation of the integral for r — > oo. 
The reason is that the pole of the integral may be very close to one of the 
branching points which is given by k and which will affect the steepest descent 
evaluation of the integral. This was the point overlooked by Sommerfeld in his 
earlier works according to JjJ, |22j though Sommerfeld mentions the problem 



44, p. 258] later but refers for further details to a paper by H. Ott 35]. In 
14 . p. 932] he also acknowledges comments from F. Noether and V. Fock 
on the perils of the approximation procedure which in the end means that 
in practice the "surface wave" cannot be separated from the space wave ("Es 
diirfte keine Bedienungen geben, unter denen sich der Oberflachenwellentypus 
P rein ausbildet und den Hauptbestandteil des Wellenkomplexes darstellt"). 
However if the ground is covered by a dielectric layer (such as ice over the 
sea) then there may indeed arise trapped surface waves in the layer as pointed 



13 They claim to have settled the issue "in the present paper by proving in a quite general 
way that this surface wave cannot be included in the said dipole radiation and by pointing 
out a thus far hidden error in Sommerfeld's computation" [22J, p. 807, abstract]. However, 
the problematics of the asymptotic evaluation of the integral involved was early recognized 
by e.g. F Noether and V Fock and communicated to Sommerfeld. 
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out by Wait [500- The branching points mentioned above arise because of 
the square root expressions when the integral is evaluated using the residue 
calculus. In this procedure the path of integration is closed by an infinite 
half circle in the upper half plane inside which there will be one pole and 
two branching points (at k and fog). Cutting up the plane along lines from k 
and k E to k + ioo and k E + ioo will make the integrand single- valued and the 
residue theorem becomes thus applicable. For further discussion of this topic 
see the references that have been listed above. Returning to the discussion 
of the Zenneck waves we may envisage that they can appear if a plane wave 
travelling along the x-axis hits a dielectric/conductor at x = 0, which extends 
to z < and x > 0. Part of the wave scatters but part of it will continue 
along x > in the half-space z > and is expected to approach the form of 
a Zenneck wave far from x — 0. 



14 The question of the presence of surface waves may be compared to similar problems 
of identifying particle states and resonances in quantum mechanics; investigations which 
too involve scrutnizing the poles in evaluating the "transfer functions" of the ip-w&ve. 
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A Some definitions and results from vector analysis 

Vectors are denoted by A = (A x , A y , A z ) = A x ~k + A y y + A z z in a Cartesian 
coordinate system; A = (A p , A v , A z ) in a cylindrical coordinate system; A = 
(A r , Aq, Ay) in a spherical coordinate system. Note that the ^-component 
appears in different order in spherical and cylindrical coordinate systems 
though it has the same geometrical meaning in both systems. 



A B = A X B X + A y B y + A Z B Z (dot product) 

A x B = (A y B z — A z B y , A Z B X — A X B Z , A x B y — A y B x ) (cross product) 
A x (B x C) = (A • C)B - (A • B)C (Lagrange identity) 

fd<P d<P d<j>\ 

v< ^= (^>~7r~>ir~ J ( _c y L coord -) 



dp p dip dz 

.' d<j) 1 d(f) 1 , , 
W = »!> I^T' —T TJTIT (-spher. coord.) 



dr' r d9' r sin # d</? 
„ . <9A„ <9A 2 

V • A = — h + -tt - (divergence) 

ot dy dz 

V • A = -- ^-^ + -- ^ + -cyl. coord. 

p dp p dip dz 

„ . 1 9(r 2 A r ) 1 d(sinM e ) 1 31 , , 

V • A = - V r; + , V ^ + , y -spher. coord. 

r z dr r sin 9 da r sin # 

\ <9z ' dz dx ' dx dy ) 

V x A = _ d At d AfL _ ^ 1 d (P A 4>) _ ^dA p \ l cQOTd 

\p dip dz ' dz dp ' p dp p dip J 

1 /9(sin04p) &4^\ 1 dA r ld{rA v ) 



V x A = 

1 rsin# ^ d9 dz ) ' r sin 9 dip r dr 

1 fd(rA v ) dA. 



r 



[-^---Qf)) (-spher. coord.) 

SP_ cP_ 

dx 2 dy 2 dz 2 



d 2 d 2 d 2 

?2 = V • V = — T + ^-r + -^-t- (Laplacian operator) 



1 <9 2 

□ = V 2 — 7717 (wave operator) 

c 2 at 2 

V x (V x A) = -V 2 A + V(V • A) 

V • (A x B) = B • (V x A) - A • (V x B) 
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(t A ds= / V x A ■ dS (Stokes-Green theorem) 

'as Js 



A ■ dS = / V • A dV (Gauss' theorem) 

av Jv 

I (FVG - GVF) • dS = [ {FV 2 G - GV 2 F) dM (Green's theorem) 
Jav Jv 



'av Jv 

„ikr 

(V 2 + k 2 ) - — = -6{t 

v 1 A-nr 



J e~ ax2 dx = ^ (»(a) > 0) (Gauss integral) 

The scalar Laplacian V 2 can be calculated for cylindrical and spherical 
coordinates by using the identity V 2 = V • (V0). For vector functions the 
Laplacian operator can be evaluated using the identity V 2 A = V (V • A) - 
V x (V x A). The Dirac delts0 5(x) is a so called functional defined by the 
property 

f(x)S(x — xo)dx = f(xo), (Dirac delta.) 



R 



for any function /. The Fourier representation of the Dirac delta can be 
expressed as 



S( x ) = — I e ikx dk. (168) 
27r Jr 

We may extend the Dirac delta to vectors by defining (in Cartesian co- 
ordinates) 8(r) = S(x)5(y)5(z). 

In the text we have often used complex notation for the EM fields. Thus, 
an harmonic E-field is written^ 

E(r,t) = ~E(r)e~ iuJt . (169) 
The physical field will then corresponds to the real part of the "phasor" 

PHD, 

E phys (r,t) =&(E(r,t)). 



15 It became well known through Paul Dirac's Principles of Quantum Mechanics (1930) 
but the concept had been introduced in electrical engineering by Oliver Heaviside (1850- 
1925) from where Dirac picked it up. 

16 We here adhere to thephysics tradition using the time factor exp(— iut), while many 
engineering texts such as [3J] assume a time dependence of the form exp(io;t) (also often 
denoting \/— 1 by j instead of i). 
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When computing dot-products and cross-products for complex fields, we 
have to use complex conjugation in order to obtain the corresponding physical 
quantities. For instance, E phys • E phys is evaluated as 

-E • E*. 

2 

The factor 1/2 comes from the fact that a time average is implied (RMS- 
value), which for real fields amounts to the factor 



- J (co S (cut)) 2 dt=-. (T = 2n/u.) 
Similarly the Poynting vector is expressed as 

S = \ 3? (E x H*) , 

and the power density V = E • J becomes likewise in terms of phasor 
quantities 



V = - ft (E • J*) 



B Tables 



Constant 


Value 


Unit 


Legend 


eo 


8.854187817 -10" 12 


As/Vm 


Permittivity of vacuum 




4tt • 1(T 7 


Vs/Am 


Permeability of vacuum 


Vo 


376.7303 


Q 


Wave-impedance of vacuum 


c 


2.99792458 -10 8 


m/s 


Velocity of light in vacuum 




1.60217733 -10~ 19 


As 


Electron charge 



Tab. 1: Constants 

Here cq and /xo are exact values. These and t]q and are connected by 

1 



c 



Vo = \ — , 

eo 



from which we obtain 



?7o = c 0/ uo ~ tt120 VL w 377 Vt. 
Sometimes tables list the "loss tangent" defined by 
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Material 


e 


a (Q my 1 


u 








1 
1 


1 act' oyhn^rn 


1 4 


10~ 7 


1 
L 




Q 

o 


in- 11 


1 
1 


uiy duck 


t± 


u.ui u.uz q.o un/j 


1 
1 




7 ^ 




1 

J. 


w nn H 
vvouu. 


Q 



in-13 in-4 
1U 1U 


1 
1 


1 ImifT QO TIT* 1 T~\ 1 T T\\jr\ An 1 

uougicto in ipiywoou.) 


1 R9 


n nziQ to ^ 


1 
1 


>ai Icon/Hi 
Dull I btlllU. J 




i n -5 


1 
1 


\ A ■ \ T / A Y ' | 1 

VVdXt!! ^bt!d ) 


ou 




1 
1 


\X7r\ Inv f 1 Q 1/A \ 

VVdXtil ^IdKt; ) 


oU 


in -3 to 'ir^'i 


1 
1 






^ -1 0~ 4 CO TtH?"! 


1 

J. 


Snow 


1.5 


10~ 3 


1 


Air 


1 





1 


Human tissue 


70 


0.2 


1 


Iron 


NA 


10 6 


5000 


Copper 


NA 


5.8 -10 7 


1 


Aluminum 


NA 


3.5 -10 7 


1 



Tab. 2: Electric properties 



e" a 



tan 5 

e' e u; 

That is, e" = e' tan 5. Electric properties depend generally on temperature 
and frequency. The complex permittivity of water, for instance, can be quite 
well represented for / < 50 GHz and 20° C by the Debye model (quoted by 

H) 

. e s - . a 

e = 6oo + ; h % , 

1 — iujt ioe 

where = 5.27, e s = 80.0, r = 10~ n s, and the final term accounts for 
the effect of salt if present (salt water having a conductivity around 3-5 
(ftm)- 1 ). 

About wood note that is usually anisotropic due to the fibers. Fields po- 
larized along the direction of the fibers pass more easily through the material 



37|. For human tissue there has been reported [41fl the following expression 



for (real) permittivity and conductivity, 
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e ' = 5 + TTTIW (170) 

a = ao+ 60A 1 + (1.5/A)*' (171) 

where V ~ 5 ("volume fraction occupied by macromolecular compo- 
nents"), cr ~ 1/70 (cm f2) _1 (observe the unit) and A is expressed in cm. 
These values are to be considered as averages since the properties vary with 
the type of organ and tissue (blood, bone, muscle, brain, etc). 



C Bessel functions 
C.l Bessel J-functions 

While cos- and sin-functions are associated with rotational symmetry in the 
plane (the circle), Bessel functions are associated with rotational (cylindrical) 
symmetry around an axis in the 3-dimensional space. Bessel functions come 
typically into play when we have to solve the Laplace equation V 2 w = in 
space for a system with cylindrical symmetry. The Laplace operator is given 
in cylindrical coordinates as, 

V 2 1 d f du\ 1 d 2 u d 2 u (172) 
U pdp\ dp ) p 2 dip 2 dz 2 

Consider the eigenvalue equation V 2 w = au in the case a = -1. Appar- 
ently u = exp(zy) = exp(ip sin tp) is a solution. We develop it as a series in 
exp(z^), 



e «p sirup 



Y,upy nv (173) 



which defines the Bessel functions J n of integer order (n = 0, ±1, ±2, • • • ). 
Inserting the expression (11721) into V 2 w + u = we find that the J n satisfies 
the equation, 

1 d ( dJ n (p)\ ( n 2 \ , s , . 

'p^t 1 + i 1 -^ )JM = o. (174) 



pdp \ dp J \ p 2 
From Eq. (11731) it follows also that J n can be defined via the integral 

f2- 



Up) = ^~ r ' e'^e-^dip, (175) 
27r Jo 
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and that we have a series expansion (for n > 0), 



„2 



= g t!r (n + t + i) (-Tj ( " >0) - (176) 

We have in Eg. (11761) employed the Gamma- functionary) = J °° s t ~ 1 e~ s ds 
> 0) which for integer values m > satisfies T(m) = (m — 1)!. Eq. (11761) 
can in fact be extended by replacing n by fractional values i/ to Bessel func- 
tions J u of fractional order which also satisfy the Bessel equation (11741) . As 
a historical note we may mention that Daniel Bernoulli studied the series 
expansion of what is now called Jo already in 1738 in connection with the 
problem of determining the shape of an oscillating chain (catena) supported 
from both en ds [52l . sec. 1.2]. Using the method of steepest descent or sta- 
tionary phase [44j, sec. 16.E] |l2|, sec. VII. 6] one can show that the asymptotic 
value of J n (p) for p — > 00 and fixed n is given by 



Jn{p) « J— cos(p-(2n + l)~) . (177) 
V Tip V 4/ 



2 

Tip 

Thus for large p the function J n (p) behaves as the cos-function. One also 
sees from this that J n (p) has an infinite number of zeros, as already pointed 
out by Daniel Bernoulli in case of J . Below we will also allow J n (z) to 
take complex valued arguments z. From small p we have from Eq. (jl76l) the 
asymptotic relation, 

Since exp(i(y) = exp(i(psmip) = J2 n Jn(Cp) exp(m^) it follows that 
J n (Cp) satisfies the equation, 

SWM'- -?)■«<»- ^ 

Especially we have in the case n — 0, 

^ + ^ + ?)m<p)=0- (180) 



17 A more general definition is (L Euler 1729) 



Y{t) = lim 



n^oo t(t + l)---(t + n) 

which is valid for all points t in the complex plane except for and the negative integers 
(the poles). 
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This implies that for any function F(() 

n= /°°F(C)Jo(Cp)e ± v / ^^C (181) 
Jo 

satisfies the wave equation V 2 n + k 2 U = 0. This fact was employed in 
section [5.4.11 in order to write exp(i£r)/r on the form (11811) . One makes 
the ansatz that there is a function F such that (I181H is satisfied in case of 
II = exp(ikr)/r. One studies first the instance z — 0, 



e ikp 



P 

Using the general orthogonality property 



poo 

/ F(()J ((p)d(. (182) 
Jo 



POO 

/ J n ((p)J n (rp)pdp = r _1 5(C - r), 
Jo 



(183) 



of the Bessel functions one can invert Eq. (11821) obtaining [44j, p. 243] 



/•OO 1 

C-^(C) = J o e^J (Cp)dp=-j= 



v^P' (184) 

Once we have the result for z = (I182p we obtain the general solution 
by inserting the factor exp(±^/£ 2 — k 2 z) as in Eq. (I148j) . 

The orthogonality relation Eq. (|183l) can be demonstrated by a clever use 
of the identity 

f( x > y } = JJ f(^ r l)$(x-0 5 (y-v)d£dri, 

inserting the representation Eq. (11681) for the Dirac deltas and then going 
over to polar coordinates assuming the special form f(r, <p) = g(r) exp(intp) 
for the function / (44, sec. 21]. That J n (rp)J n ((p)pdp = for r ^ ( 
can be seen also directly by multiplying the Bessel differential equation for 
J n (jp) by J n {(p) and vice versa and subtracting the expressions. 



C.2 Hankel functions 

As we found above the Bessel function J n resembles the cos-function. It is 
often more useful to employ the complex exp-function than the cos-function, 
for instance when dealing with time-harmonic fields. It likewise exists a 
Bessel counterpart to the exp-function that may be preferable to the J- 
functions in some circumstances. These are the so called Hankel functions 
Hu and . In order to arrive at these one may start by generalizing the 
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integral representation fllT6f) to one that covers also the case of non-integral 
orders (Schafli 1871), 



Up) 



1 

2^ 



ipcos z+iv(z—n 



' 2 >d 



Z. 



Ho 



The integration is along an infinite path 
Wq in the complex z-plane, chosen such as to 
make the integral convergent for real positive 
p. The standard choice for W is the path 
-f + zoo ->• -f + iO -> ? f + zO -> f + zoo. 
In case of an integer order v = n the expres- 
sion (11851) reduces to the earlier one (11761) . since 
the integrals along the vertical sections of the 
path Wq cancel each other in this case. In the 
adjoining figure we introduced two other infi- 
nite paths W\ and W 2 which also define two 
convergent integrals (introduced by Sommer- 
feld in 1896) and the corresponding functions 
called the first and second Hankel functions, 



Wo 



(185) 



3tt 
2 



Wo 



71 



71 



Hi 



W 2 



Jpcos z+iu(z-w/2) 
Jp cos z+iu(z-n /2) ^ 



(186) 



Now we have formally Wo = W1+W2 since the integrals along the negative 
imaginary axis cancel each other. We therefore obtain, 



Up) = \{h£\p) + h?\ p )) 



(187) 



This can be compared to cos<^ = (exp(iip) + exp(—i<p))/2. Indeed the 
asymptotic form as p — > 00 is given by 




3 iO-(l/+l/2)7T/2) 



-i(p-(f+l/2)7r/2) 



(188) 



For p — > we have instead 
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4%) 



■ 2, 
-%— in p, 

7T 



(189) 



TM 



7T 



(9fo/ > 0). 



One property following from the definitions (I188j) is that for real p and v 



Hi l \ P ) 



(190) 



where the bar denotes complex conjugation. The Bessel function corre- 
sponding to the sin-function is called the Neumann function and is defined 

by 



(191) 



The Hankel functions play an important role in Sommerfeld's theory of 
the dipole over ground and e.g. in the treatment of EM-fields along circu- 
lar conductors. Using the decomposition A187B of the Bessel function J in 
Eq. (1157p one can transform the /^-integral into an /^-integral. 

This makes it possible to apply the 
calculus of residues and evaluate part of 
the integral in terms of a pole of the inte- 
grand. An important step is the relation W 2 



H®{-p) = -H™{p) 



(192) 



H$\-p) = -HX\p) 



(2), 



which is a special case of the Um- 
laufrelationen for the Hankel functions 
44 p. 314]. Note that the integration 




paths used in order to define Bessel func- 
tions need to stay in the limits of \z\ -^00 
in the checker board patterns shown in 
the adjoining figure in order to ensure convergence for p 
the checker pattern is shifted by n along the real axis), 
though be deformed askong as they do not cross borders and still define the 
same functions. Thus W2 is an allowed version of W2 described earlier. Now 
H o\~p) (P > °) is defined using the path W\ translated by 7r along the 
real axis, which we denote W\ + n. Since — pcosz = pcos(z ± 7r) we get 



> 0; for p < 
The paths may 
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-^o \~P) = (V 71 ") J w exp(ip cos(z))dz where W\ + tt is shifted by tt for z > 
and by — n for z < 0. This gives a path equivalent to W<i except that it is 
traversed in the opposite direction. 

D The cable equation and impedance 

Waves do not travel only in open space, but also along cables. Thus, for 
high frequency currents, circuit characteristics can no longer be treated ac- 
cording to the usual point-to-point models. For instance, the lengths of the 
connections may affect the circuit properties. One central property of circuit 
elements, cables and loads, is that of impedance Z. Generally speaking, if 
we feed a voltage V — V • exp(-iut) into a load and measure a resulting 
current I — Iq • exp(-iujt), then the impedance (at the given frequency) of 
the load is defined by the quotient 




We will consider a discrete model of a cable made up of a series of induc- 
tors L and capacitors C (in parallel). 



Fig. 4: Discrete model of a transmission line as a series of inductors L and 
capacitors C. 

These elements are defined by their impedances Zl = —iuL and Zc = 
1/iuC; that is, the voltage over an inductor is given by V = L dl/dt and 
over a capacitor by V = Q/C. If one applies the laws of Kirchoff and Ohm 
to the currents and voltages over the n — 1, n, and n + l:th elements one 
obtains the equations 




L 



-CV r 



—LI n — V n+ i — V n . 



n 



n 



I, 



n—li 



(193) 
(194) 
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Letting the spacing Ax between the units go to zero we may replace the 
difference by derivatives, 



~CV=f x , (195) 
dV 

-CI = w (196) 

where C = C j Ax is the capacitance per unit length, and £ = Lj Ax the 
inductance per unit length of the cable. This leads to the "cable equation" 

d 2 V 1 d 2 V , 1 , . 

with u = -=. (197) 



dx 2 u 2 dt 2 " ' ' y/CC 

This is a wave equation for EM-waves propagating along the cable with 
the velocity u. Inserting solutions of the form V = V ■ exp(i(kx — cut)), I = 
J • exp(i(kx — ui)) into (11951) . where k = u/u, we obtain for the impedance, 

Thus, the quantity Z = ^J~^ gives the cable impedance for this model. 
One can also add resistive elements to the model, but the general features 
are already apparent in this simple model. Standard coaxial cables have 
impedances of 50 Q or 75 f2, but one may note that these values apply only 
in some restricted frequency range. 

In the above case we considered an infinite long cable. Suppose instead 
we have a cable with impedance Z of length / which is ended by a load Z L . 
What will then be the impedance of this system? Assume the feeding point 
is at x = 0, and the load is at x — I. If we feed a voltage Vq ■ exp(— iut) at 
the point x = 0, then the solution for the system is of the form 



V(x, t) = ae i{ - kx - ut) + fte^-**-"*) , (199) 
I(x,t) = ce i < fcB - wt) + de i{ - kx ~ ut) . 

The terms with — ikx in the exponents correspond to the reflected part 
of the wave. The impedance will thus be given by the quotient 

_ V(0,t) _ a + b 
~ 7(0, t) ~ c + d' 
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The coefficients a, b, c, d are determined by 
the condition V(l,t) = Z L ■ 1(1, t), and using 
(11951) (which yields the relations a = Z c, b = 
—Z d). After some algebra one obtains finally 



Z 



Zl — iZ tan(fcZ) 



(200) 



' Zq — iZl tan(fcZ) 

We have the interesting result that the 
impedance depends on the cable length I and 
the wavelength A through the term tan(27r//A). 
For instance, for a quarter wave cable, I = A/4, 
short-circuited case (Zl = 0) we obtain 



we get Z = Zq/Zl- In the 



Z = —iZ taa(kl). (Short circuited case.) 

In this case we have infinite impedance for quarter wave cables. This 
property is used in microwave chokers which are quarter wave cavities de- 
signed to trap unwanted leakage of EM-waves e.g. through the door slits 
microwave ovens. 

For an open ended cable (Zl = oo) we have instead 

Z = iZ cot(kl), (Open ended case.) 

and this corresponds to an infinite impedance for a half wave cable. We 
can calculate the rms power dissipated by the load Zl from (assuming Z is 
real) 



^ = ^(0n0) = 2^' , - ,a 



V 2 



(H 2 -|6| s 









2 


1 - 


z L - 


-z 






l-Z 





2Z 



Zl-Zq j2M 
Z L +Z C 



2 ' 



The quotient \b/a\ describes the fraction that is reflected, and it is ob- 
tained from 



- = jr ■ e 1 ™- (201) 
a Z L + Z 

From this we infer that there is no reflexion when we have impedance 
matching, Zl = Z . In this case the energy transferred to the load is sim- 
ply V 2 /2Z . The term impedance matching is also used for point-to-point 
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circuits where reflexion plays no part. If we have resistances Rq and Rl in 
series and apply the voltage V, then the power dissipated by the resistance 
Rl becomes 



Rr. ■ I 2 



R 



L • 



V 



v Ro + Rl , 

which for fixed V and Rq attains its maximum when Rl = Ro- Thus, for 
this case, in order to deliver maximum of power to the load, its impedance 
must be matched to the circuit impedance Ro. If we repeat this exercise 
with the complex impedances Z and Zl we obtain the matching condition 
Zl = Zq. These results can be applied to the case where we have an RF- 
circuit with impedance Z coupled to an antenna with the impedance Zl- 
From (12011) we get for the magnitude of the reflexion coefficient, 

IPI = f^! 2 ■ (202) 

Two antenna design parameters are defined in terms of r: the Voltage 
Standing Wave Ratio, 



and the Return Loss, 



VSWR 



i + H' 



(203) 



ft = 20 log(|p|). 

Designing the antenna conditioning one has 
to make a choice between minimizing the reflec- 
tion (Z = Zl) and maximizing the radiation 
power (Z L = Zq). One wants to get rid of the 
reflexions because they can interfere with RF- 
circuit operations. Above we considered a dis- 
crete model of a transmission line. The straight 
coaxial cable can be easily solved too because 
of the cylindrical symmetry. Thus assume the 
center wire has the radius a and outer conduc- 
tor an inner radius b. The transversal electric 
and magnetic mode (TEM) correspond to the 
case where the electrical field E = (E r , 0, 0) 
is radial while the magnetic field is azimutal, 
H = (0,^,0). The Maxwell equations of in- 
terest in the region a < r < b become, 



(204) 
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dE r „ , N 

-r£ = iooB 9 , (205) 

9 - iuB^ (206) 



dz 

■n \ 

0. (207) 



1 d{rH v 



From this follows that 



r dr 



d 2 E r , 9 
+ A; E r 



where k = u^/eeoJIJIo- It follows from the third equation that H v oc 1/r. 
More precisely we have ^ H ■ dr = I where the integral is along a circle of 
radius a < r < b encircling the inner wire which carries a current /. Thus 



which in combination with the equation (12051) gives 



////To I I 
E r = J — — = v: 



The potential V between between the inner and outer conductor becomes 



therefore 

rb 



ril fa 
E r dr = — m — 
2vr \b 



Finally the characteristic impedance Z — V/I of the cable becomes 

2vr 



Z=^ln^ 



E Chipcon WCR2400 



The Chipcon device CC2420 [10] operates in the 2.4 frequency band (A = 
12.5 cm). IEEE 802.15.4 defines 16 channels in 5 MHz with the frequencies 
f k = 2405 + 5 • (k - 11), with k = 11, 12, . . . , 26. The effective data rate is 
250 kbps (2 MChips s -1 ) and the device uses a coding with 4 bit symbols in 
32 chip spread sequences. The RSSI values are determined from the average 
over 8 symbols corresponding to a time interval At = 128 //s. The signal 
strength is determined from the automatic gain control factor in the signal 
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amplifier part (variable gain amplifier, VAG). From the bit rate we gather 
that every bit has a 4 /xs time window. Since the velocity of light c is about 3 
x 10 8 ms" 1 we get c x 4 /is ~ 1200 m for the distance that the radiation can 
travel during the time window which means that the reflected and refracted 
parts (from objects within 600 m which by far exceeds the normal range 
of the devices) can contribute to the bit-signal, either by constructive or 
destructive interference. The corresponding chip time window is 1 /js which 
translates into a 300 m distance. In our measurements we used a standard 
coaxial A/2-dipole antennas model WCR-2400-SMA. The dipole antennas 
approximate well a circular radiation pattern in the horizontal plane @. 




Fig. 5: Variation of RSSI (radius corresponds to absolute RSSI value) with 
angular orientation of the receiver with shielding (circle) and no 
shielding (filled circle) in outdoor setting. The receiver-transmitter 
separation was 10 m. Origin of the diagram corresponds to RSSI = 
-50, and the radial intervals to 10 RSSI units. 

However, the antenna + RF device results in orientation effects appar- 
ently because of the PCB affects the radiation pattern of the antenna in 
the horizontal plane. Shielding the RF device using an appropriate metal- 
lic enclosure ("Faraday cage") reduces the orientation effect significantly as 
shown by FigJHl The orientation effect is similar in the transmitting and 
receiving mode. One may also note that we get a better RSSI value in the 
shielded case. Cylindrical enclosures of diameter 90 mm and height 40 mm 
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Fig. 6: Shielding box. 



were custom made in 4 mm aluminum. The enclosure housed the PCB plus 
the battery, with the antenna fitted to the center of the top surface (top and 
bottom plates 1 mm thick). 

The enclosure + antenna naturally will be expected to have a different 
input impedance than the antenna alone. However, there seems not be any 
readily available treatments of how such a cylindrical enclosure modifies the 
impedance. This problem as well as impedance measurements using a net- 
work analyzer will be left to a separate study. 



F Reciprocity theorem 



A basic property used in antenna measurement [24| is the so called reciprocity 
theorem which says that an antenna has the same field pattern both in the 
receiving and transmitting mode. Therefore one needs to test it only as either 
receiver (usually the simpler alternative) or as a transmitter. More generally, 
suppose an antenna Ti at rj is transmitting radiation at the frequency / and 
is intercepted at r 2 by an antenna T 2 , then field strength measured with T 2 
is the same as would be measured Ti if T 2 is transmitting with the same 
energy and at the same frequency /. If we imagine the antennas as part of 
of circuit [s| with voltages and currents related by 



V 1 = Z n h + Z 12 I 2 (208) 
V 2 = Z21I1 + Z22I2 

then reciprocity boils down to the symmetry of the mutual impedance, 
Z12 = Z21: A given current in one antenna causes the same (open circuit) 
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voltage in the other antenna irrespective which of them acts as the transmit- 
ter or receiver. 




O — I I — 1 — I I — o 



J V2 



Fig. 7: Antenna system as a two-port network 



The electromagnetic reciprocity theorem goes back to the work by H A 
Lorentz in 1895-1896 and an account is given by Sommerfeld in p. 653- 
663] (which however rests on the assumption that the antenna length L is 
insignificant compared to the wavelength A). It is built on the fundamental 
equation 

V • Ei x H 2 - V • E 2 x Hx = E 2 • Ji - Ex • J 2 (209) 

which is obtained from Maxwell's equations for (E 1; H 1; J x ) and (E 2 , H 2 , J 2 ) 
in the harmonic case (time derivatives of the fields are replaced by multiplica- 
tion by — ioo) assuming we have a linear and isotropic medium (D*. = eeoE^; 
B fc = /i/i H fc ; k = 1, 2), 



V x E fc = iuB k , (210) 

V x Hfc = — iuDk + Jfc- 

We may take the index 1 to refer to the case when device 1 acts as sender, 
and the index 2 to refer to the case when device 2 acts as sender, everything 
else remaining similar. The fields have to satisfy the boundary conditions on 
the conductor surfaces which is why we have in general Z 12 7^ 0@ If one 
integrates Eq. (12091) over a spherical volume containing the antennas and the 
devices we get for the left hand side, by using Stokes' theorem, the surface 
integral 

18 If the indexes 1 and 2 refer to different instances, then E& will be the total electric 
field and the boundary condition has to be applied to it. In this case Jfc will differ from 
zero not just for the sending antenna but also for the receiving antenna because of the 
induced current. If the indexes refer to the same instance but where Ji is the current of 
antenna 1 alone (and vice versa for J2) then the boundary condition has to be applied to 
the total field E x +E 2 . 
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ExxHs-rfS-/ E 2 xHi-dS (211) 

which goes to zero as the radius r of the sphere goes to infinity. The 
reason for this is that the field at a far away distance approaches locally a 
planar electromagnetic radiation for which H = ±k x E/ (kr/). Inserting this 
in Eq. (121 II) we see that the difference vanishes. Thus we get from Eq. (l211j) 

E 2 • J 1( 1V = J Ei • 3 2 dV. (212) 

If the antenna 1 acts as a transmitter it emits a power Pt = l/2-9?(Vi I* ) = 
1/2 • 3?(Zn|/i| 2 ) from which we deduce that Z\\ = Ru + iX\ where R\ s is 
the radiation resistance of the antenna while X\ is called its reactance, and 
mutatis mutandis for Z?o. If the receiving antenna 2 is furnished with a load 
Zl the current I 2 will bd_j —Z 2 \I\j [Zt + Z 22 ) and the received power hence 
given by Pr = 1/2 ■ \Z 2 xIi/(Zi + Z 22 )\ 2 ?R,(Zl). This attains the maximum 
for a matched load Zl = Z% 2 in case which we obtain the ratio of received to 
transmitted power as, 



Pr \Z\2\ 



12 



Pt AR\ s R 2s 

which is, as seen, symmetrical in the in the indexes 1 and 2 if Z\ 2 = Z 2 \. 

We will consider the reciprocity relation in a bit more detail using as 
surfaces of integration the surface S\ which follows the surface of the antenna 
1+device 1 except where it crosses a section of the the coaxial cable (see the 
right part of Fig. (ED and and the surface S 2 defined mutatis mutandis 
for the second antenna + device. One thus integrates over a volume V with 
the boundary Si + S 2 which excludes the interior of the antennas+devices 
and where we therefore have Ji = 3 2 = 0. One obtains then the reciprocity 
theorem in the form 



E 1 x H 2 ■ dS - I E 2 x Hi • dS = 0. (213) 



<S 1 +S 2 JS!+S 2 

Here, for instance, the integral J s Ei x H 2 • dS becomes J A E x x H 2 • dS 
where A 2 is a cross section of the coaxial cable interfacing the antenna 2, 



19 Indeed, the antenna resistance has to be included in the total impedance because the 
antenna reradiates also in the receiving mode. The present analysis assumes that the 
distance between the antennas is large compared with the wavelength A and thus Z12 
small in comparison with Zu and Z22', that is, the coupling between the antennas is weak 
- see further 0, sec. 4.2]. 
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because elsewhere the electric field is orthogonal to the surface meaning that 
Ei x H 2 • dS = — H 2 ■ E x x dS = . In this integral H 2 is the field in the coaxial 
cable generated by the impressed current J 2 while Ei refers to the field in 
the coaxial cable induced by the impressed current I\ in the antenna 1 via 
the transmitted field impinging on antenna 2. In the integration is only 
over S% and the volume V thus contains the antenna 1 and we will therefore 
have a current term J v E 2 • JidV on the right hand side of Eq. (12131) . In 0, 
Eq.(15) and the preceding one] I\ apparently refers to the current induced 
by the field Ei in antenna 2 (a point left somewhat unclear in 




Fig. 8: Sketch of the dipole antenna link to the coaxial cable. On the right 
the arrows indicate the flow of field energy in the receiving mode. The 
dipole conductors steers the flow along the conductor surface into the 
coaxial cable. In the transmitting mode the energy flow is outwards 
from the coaxial cable and directed to the surrounding space by the 
antenna surface. On the left there is a schematic outline of the surface 
S of integration. The surface S follows the dipole antenna and the 
receiver/sender along the conducting surfaces where the tangential 
fields are zero, except where it forms the cross section A. 

For a TEM wave in a coaxial cable we have the solution 
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H v {r) = ±— , (214) 
E r (r) = r] 1 



2nr ' 



where I is the current of the coaxial cable. The electric field is thus radial 
and orthogonal to the magnetic field which is along the angular direction. 
Inserting this solution into f A Ei x H 2 • dS we obtain 



Ej x H 2 ■ dS = -^-I^ 2 I 2 [ -2-nrdr 
(2tt) 2 J a r 2 



2/T \a J 

Here a is radius of the inner wire of the coaxial cable and b the radius 
of the outer conductor, while Z c defines the impedance of the cable. We 
use 7x->2 to denote the current induced in antenna 2 by the field radiated by 
antenna 1 with the impressed current I\. The evaluation of f A E 2 x H 1 ■ dS 
will give the same result except being of the opposite sign (this has to do 
with the orientation of the fields and the fact that the impressed and induced 
fields go in opposite directions in the coaxial cable). This leads finally to the 
relation 

2Z c J 1 _> 2 / 2 = I E 2 -J 1 dV. 
Jv 

If the receiver circuit is coupled to a matched load Rl = Z c then 2Z c /i^ 2 
will correspond to the induced voltage V 2 = Z21I1. From Eq. (12121) and re- 
peating the above considerations for antenna 1 we obtain Z 12 = Z 2 \. The 
induced voltage in antenna 2 can, according to these results, be written aS 

V 2 = ^ / E x ■ J 2 dV. (215) 
h Jv 



G Measurements of the reflexions from the ground 

The measurement procedures were designed such as to obtain knowledge 
how the RSSI depends on the distance the receiver and transmitter, and on 

20 When Ei refers to the total field then the result corresponds to closed circuit voltage, 
when Ei refers only to the incoming field then the result corresponds to the open circuit 
voltage. 
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the environment. The devices were typically attached to wooden supports 
(poles) of height 1 m. Power sources were 9V batteries. The transmitter 
was configured to send 20 packets per transmission. The transmissions were 
intercepted using a "sniffer" device (Chipcon Packet Sniffer) connected to a 
laptop computer. The data was saved on the computer for later extraction 
of the RSSI values, packet information and other data. For a majority of 
the tests the channel nr 20 was used and the power set to level 11 (OxAOB) 
corresponding to -11 dBm (see Tab. (jHJ)). When calculating the average RSSI 
value the maximum and minimum values were dropped in order to eliminate 
possible outliers. Alternatively we used the median value. The RF-circuit 
was shielded by putting the whole PCB plus battery into a metallic can of 
90 mm diameter and 40 mm height. 



PA level 


Output power (dBm) 


Register 


31 





OxAOFF 


27 


-1 


OxAOFB 


23 


-3 


0xA0F7 


19 


-5 


0xA0F3 


15 


-7 


OxAOEF 


11 


-10 


OxAOEB 


7 


-15 


0xA0E7 


3 


-25 


0xA0E3 



Tab. 3: Power level assignments 



In case of a smooth ground the two-ray (mirror) model can be fitted 
quite well to the data, see FigH In the figure we compare measured data 
(RSSI values) with data (power in dB) computed using the two-ray model 
described above with the dielectric permittivity e set to 3 providing a good 
fit. Measurement points were sampled more densely (Ar = 10 cm) where we 
expected the changes to be largest. In this case the Brewster angle 9b will 
be 30°. The Brewster distance is about 3.9 m, and the breaking point is at 
r « 42 m. 
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Fig. 9: RSSI and power dB vs distance, measured data (circle) and data 
computed using the two ray model (+) with e = 3. The tops of the 
antennas were about 1.14 m above the ground (sand). At the time of 
measurement the temperature was -5°C and the ground was covered 
with a ca 5 cm layer of fresh snow. 



hi(m) 


Correlation 


Error (stdev) 


Offset (mean, n) 


0.64 


0.99 


0.95 


49.4 


1.14 


0.99 


1.15 


48.3 


1.64 


0.98 


1.30 


49.3 



Tab. 4: Statistics of model-data comparison 



Similar measurements were made with h\ = 0.64 m, 1.64 m, while hi 
remained 1.14 m, with equally good fits between model and data. Tab.(Jlj) 
shows the statistics of fitting the measured RSSI values to 10 ■ log(P) — /i, 
where P is the power calculated from the model, and \i is the "offset" value 
between measured RSSI and model power (in dB). The difference (error) 
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between the model and the data are seen to correspond to about 1 RSSI 
unit when measured in terms of the standard deviation. Note that the data 
points are not evenly distributed but were chosen such as to best cover the 
places where the changes in the RSSI were expected to be largest. 

Measurements were also made on a frozen river with a smooth ice of 
thickness around 10 cm. The results were modeled using the dielectric layer 
model treated in section [3751 see Fig. (fT0lh A characteristic difference when 
compared with Fig. ((9]) are the much deeper troughs of the interference pat- 
tern in case of sandwiched structure air-ice-water. In the model this feature 
is sensitive to the thickness of the ice. 




Fig. 10: RSSI and power dB vs distance measured on a frozen river. Mea- 
sured data (+) and data computed using the two ray model (contin- 
uous) with two reflecting interfaces, with e = 3 for ice (10 cm thick) 
and e = 80 for water. The tops of the antennas were about 1.14 m 
and 1.16 m above the ice. 
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H C-code for extracting RSSI values 

Below is an example of the C-implementation of the function extract () 
which has been used to extract the RSSI values from the psd-files obtained 
with the Packet Sniffer. The C-code will depend on what sort of information 
the devices are programmed to relay besides RSSI (e.g. battery power level, 
sensor data). 



#include <stdio.h> 
^include <stdlib.h> 
^include <ansi_c.h> 

/* Structure used */ 

struct sRSSIval { 

int n; // size of RSSI— array 

int *RSSI; // array of RSSI— values 

}; 

typedef struct sRSSIval tRSSIval; 

/* extract defines a function which opens 

a binary psd—file and returns the RSSI values */ 

tRSSIval extract (char *fname) 
{ 

int i , j = 0; 
unsigned char temp ; 
unsigned char dummy [121]; 
int x=0, y = 0; 
int temppu ; 
tRSSIval data , err ; 

e r r . n = ; 
err .RSSI = 0; 

fp = fopen(fname, "r+b"); 
data . n = ; 
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if (fp = NULL) { 

e r r . n = — 1 ; 
return err ; 

} 

f r ead (&temp , sizeof ( char ) , 1, fp ) ; 

x=temp ; 

f r ead (&temp , sizeof ( char ) , 1, fp ) ; 

y=temp ; 

y=y«8; 

y=y+x ; //number of packets 

if (y < i){ 

e r r . n = — 2 ; 
return err ; 

} 

data. RSSI = (int*) calloc(y, sizeof ( int )) ; 

for (i = 0; i < y; i++) 
{ 

f read (dummy, sizeof ( char ) , 26, fp ) ; 
f r ead (&temp , sizeof ( char ) , 1, fp ) ; 
temppu=temp ; 
if (temppu! = 0) { 

if ((temppu&0x40)>0) { 

temppu=temppu&0x7f ; 

temppu=temppu "0x7f ; 

temppu=temppu + 1 ; 

temppu=— 45— (temppu ) ; 

data . RSSI [ j ] = temppu ; 

j++; 

// increments true 
// number of RSSI vals 



} 

else { 



temppu=temppu&0x7f ; 
temppu=— 45+(temppu ) ; 
data . RSSI [ j ] = temppu ; 

j++; 

// increments true 

// number of RSSI vals 
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} 

} 

f read (dummy, sizeof ( char ) , 105, fp ) ; 

} 

f close ( fp ) ; 
data . n = j ; 
return data ; 
extract end 
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